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1 Introduction 

1.1 Description of the problems 

Simple transport processes, such as those observed near a spatially localized ther- 
mal contact or tunnelling junction between two macroscopically extended metals 
at different temperatures and chemical potentials, have been studied experimen- 
tally and theoretically for a long time; see e.g. [Ma]. One is interested, for ex- 
ample, in measuring or predicting energy and charge transport through a contact 
between two metals, as well as the rate of entropy production. The natural theoret- 
ical description of such processes is provided by quantum theory, more precisely 
by non-equilibrium quantum statistical mechanics. The results of experiments 
or theoretical calculations can, however, often be expressed in the language of 
thermodynamics. In this paper we attempt to study such transport processes in a 
mathematically precise way, extending or complementing results in [DFG, EPR, 
JP1, JP2, Rul,Ru2]. 

In recent years, interest in transport processes has been driven by various ex- 
perimental and theoretical developments in mesoscopic physics and the discovery 
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of rather unexpected phenomena. Among them we mention dissipation-free trans- 
port in incompressible Hall fluids [La, TKNN, BES, ASS, FST], or in ballistic 
quantum wires [vW, Bee, ACF, FP1]. In such systems, "transport in thermal equi- 
librium" and the quantization of conductances are observed. Further interesting 
transport processes are electron tunnelling into an edge of a Hall fluid [CPW, 
CWCPW, LS, LSH], and tunnelling processes between two different quantum 
Hall edges through a constriction leading to measurements of fractional electric 
charges of quasi-particles (see, e.g., [SGJE], and [FP2] for theoretical consider- 
ations). At present, these processes are only partially understood theoretically. 
Other examples are Josephson junctions and Andreev scattering [SRI, SR2], or 
energy transport in chains (see, e.g., [Af] and references given there). 

In this paper, the main emphasis is put, however, on conceptual aspects of the 
theory of simple dissipative transport processes between two quantum-mechanical 
reservoirs and on an illustration of the general theory in a simple example, namely 
transport of energy and charge between two metals, described as non-interacting 
electron liquids, at different temperatures and chemical potentials. Of particu- 
lar interest to us are connections between theoretical descriptions based on non- 
equilibrium quantum statistical mechanics, on the one hand, and on thermody- 
namics, on the other hand. Our quantum-mechanical description involves equi- 
librium and non-equilibrium states of macroscopic reservoirs with many degrees 
of freedom. We show that, on intermediate time scales, tunnelling processes can 
be described in terms of non-equilibrium stationary states (NESS), examples of 
which have recently been studied in [DFG, EPR, JP1, JP2, Rul, Ru2, BLR]. Our 
construction of non-equilibrium stationary states is based on methods of algebraic 
scattering theory and is inspired by ideas in [He, Rob, BR, BM, Ha]. Links be- 
tween quantum statistical mechanics and thermodynamics are constructed by pro- 
viding precise definitions of energy and particle currents and of entropy produc- 
tion and by deriving a suitable form of the first and second law of thermodynamics. 
The general theory developed in Sect. 2 reviews ideas and results scattered over 
numerous articles and books and represents an attempt to provide a somewhat 
novel and, we believe, rather clear synthesis. A more complete version, including 
the treatment of systems with time-dependent Hamiltonians, appears in [FMSU] 
and in a forthcoming paper. It is illustrated on the example of two reservoirs 
of non-interacting electrons coupled through local many-body interactions (Sects 
3 through 5). Examples of non-equilibrium stationary states supporting particle 
and/or energy currents are constructed with the help of a convergent perturbation 
(Dyson) series in the many-body interaction terms. The currents and the entropy 
production rate are calculated quite explicitly to leading order. This enables us to 
show that, under natural hypotheses, they are strictly positive. Onsager reciprocity 
relations are established to lowest non-trivial order in the many-body interaction 
terms. Positivity of the entropy production rate has also been established recently 
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in [AP, MO] for XY chains, and in [CNP] for wave turbulence. 

1.2 Contents of paper 

In Section 2.1, quantum-mechanical reservoirs are introduced, whose time evo- 
lution is given in terms of a one-parameter group of * automorphisms of a kine- 
matical C*-algebra of operators. Conservation laws of reservoirs are described by 
commuting conserved charges. The equilibrium states of such reservoirs are in- 
troduced and parametrized by temperature and chemical potentials. Two general 
assumptions, (Al) and (A2), are formulated. They state that the thermodynamic 
limits of the time evolution and of the gauge transformations of operators in the 
kinematical algebra of an infinite system and of the equilibrium states exist. 

In Section 2.2 we study two interacting reservoirs at different temperatures and 
chemical potentials. Each reservoir is required to satisfy the assumptions formu- 
lated in Sect. 2.1. A class of many -body interactions coupling the two reservoirs 
is introduced. It is assumed that the thermodynamic limit of the interactions and 
of the corresponding time evolutions exist. Energy and charge currents for finite 
and infinite systems are then defined. 

Connections between quantum statistical mechanics and thermodynamics are 
elucidated in Section 2.3. The entropy production rate is defined and expressed in 
terms of the currents and of thermodynamic parameters. An inequality expressing 
the positivity of relative entropy is shown to imply that the total entropy produc- 
tion is non-negative (see also [Ru2, JP1]). 

In Section 2.4, non-equilibrium stationary states for coupled reservoirs are 
introduced. They can be construced with the help of scattering (M0ller) endo- 
morphisms of the kinematical algebra of the infinite coupled system. A precise 
assumption concerning the existence of scattering endomorphisms is formulated. 
Our approach has its roots in Hepp's work on the Kondo problem [He] and Robin- 
son's analysis of return to equilibrium in the XY spin chain [Rob]. Robinson's 
ideas have been put into a general context in [BR, Ha]. The scattering approach is 
the starting point for numerous heuristic studies of thermal contacts and tunnelling 
junctions (see, e.g., [Ma]). The first mathematically rigorous implementation of 
this approach in a study of energy transport in the XY spin chain and of tunnelling 
between free-fermion reservoirs appeared in [DFG] . 

In Section 2.5, long-time stability properties of equilibrium and non-equilibrium 
stationary states against perturbations of the initial state of the coupled system are 
studied, and conditions for the existence of temperature or density profiles in non- 
equilibrium stationary states are identified. 

The general theory of Section 2 is illustrated in Sections 3,4 and 5 on the ex- 
ample of two coupled free-electron reservoirs. 

In Section 3, the quantum theory of finite and infinite reservoirs of free elec- 
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trons is briefly recalled, and a class of local many-body interactions between two 
such reservoirs satisfying the general assumptions formulated in Section 2 is in- 
troduced. 

Our main technical result, the existence of scattering (M0ller) endomorphisms, 
defined on an appropriate kinematical C* -algebra describing two infinite free- 
electron reservoirs, is established in Section 4. A similar result has previously 
been proven in [BM]. We show that, under appropriate smallness and regularity 
assumptions on the many-body interaction terms, the scattering endomorphisms 
are given by a (norm-) convergent Dyson series. As a consequence, non-equilibrium 
stationary states can be constructed with the help of a convergent perturbation ex- 
pansion. 

The results of Section 4 are used in Section 5 to derive explicit expressions for 
the energy and charged-particle currents to leading order in the many-body inter- 
action terms. These expresssions, along with the convergence of the Dyson series, 
prove that, for small coupling constants, the entropy production rate is strictly 
positive, Ohm's law holds to leading order in the voltage drop between the reser- 
voirs, with a resistance whose temperature dependence can be determined, and 
the Onsager reciprocity relations hold to leading order. 

We conclude this introduction with explicit formulae for the leading contri- 
butions to the particle current, J, and to the energy current, V, between two 
reservoirs, I and 77, of free electrons coupled to each other by a quadratic local 
interaction term with a form factor w((—k, II), (I, I)) and a brief discussion of 
the qualitative implications of these formulae. These currents are given by 

J ~ 2n [ dkdl5(\k\ 2 -\l\ 2 )\w((-k,II),(l,I))\ 2 (p n (k)- Pl (k)), 

V ~ 2tt/ dkdl \k\ 2 5(\k\ 2 - \l\ 2 ) \w({-k, II), (/, J))| 2 (p n (k) - pj(Jfe)) , (*) 

where p r (k) is the Fermi distribution of the free electron gas (r = /, II labels the 
reservoirs), and w((—k, II), (I, I)) is the interaction kernel describing scattering 
of a particle in an initial state with energy \l\ 2 localized in reservoir I to a final 
state with energy \k\ 2 localized in reservoir 77. 

If both reservoirs have the same chemical potential and the temperatures sat- 
isfy T 1 < T H then J and V are positive; particles and energy flow from the hotter 
to the colder reservoir. Similarly, if the reservoirs have the same temperature, then 
particles and energy flow from the reservoir with the higher chemical potential to 
the other one. 

Formulae (*) prove that the leading contribution to the entropy production rate 
is strictly positive, unless both reservoirs are at the same temperature and at the 
same chemical potential. 

Another consequence of (*) is that, to leading order in the interaction, and for a 
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small voltage drop, A/i = p} 1 — pj (at a fixed temperature, T, for both reservoirs), 
Ohm's law is valid, i.e., the voltage drop is proportional to the current, 

A/i ~ /?(//, T)J. 

Our calculations show that the resistance Rip 1 , T) grows linearly in T, for large 
T, it has a positive value at T = and may increase or decrease in T, at small 
temperatures, depending on properties of the interaction kernel w modelling the 
junction between the two reservoirs. 

This paper is dedicated to Klaus Hepp and David Ruelle on the occasion of 
their retirement from active duty, but not from scientific activity. Some of their 
work plays a significant role in the analysis presented in this paper. J. F. is deeply 
grateful to them for their generous support and for everything they have con- 
tributed to making his professional life at ETH and at I.H.E.S. so pleasant. 

Acknowledgements. J.F. thanks S. Dirren and G. M. Graf for useful discus- 
sions during an early stage of his efforts. We are grateful to R. Fernandez, B. 
Nachtergaele and I. M. Sigal for numerous helpful discussions. We have greatly 
benefitted from studying the works of V. Jaksic and C.-A. Pillet and of D. Ruelle 
on related problems. D.U. acknowledges the hospitality of ETH Zurich, where 
most of this work has been carried out. 

2 Elements of a general theory of junctions and of 
non-equilibrium stationary states 

2.1 Quantum theory of reservoirs 

We start our general analysis by describing "quantum-mechanical reservoirs". 
A reservoir is a quantum system with very many degrees of freedom, e.g., an 
electron liquid in a normal metal, a superconductor, a gas of atoms, or a large 
array of coupled, localized spins, but with a small number of observable physical 
quantities. It is confined to a macroscopically large, but compact subset, A, of 
physical space IR 3 . Its pure states correspond, as usual, to unit rays in a separable 
Hilbert space, H A , and its dynamics is generated by a selfadjoint Hamiltonian, 
H A , acting on the space H A . The kinematics of the reservoir is encoded into an 
algebra, A A , of operators contained in (or equal to) the algebra of all bounded 
operators on H A l . The time evolution of an operator a on H A is given, in the 

'The algebra A A is sometimes called algebra of "observables", a commonly used, but unfor- 
tunate expression. 
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Heisenberg picture, by 

af(a) := e ftffA /* a e -i«* A /» (2.1) 



and it is assumed that a A (a) G *A A , for every a G *A A . 

There may exist a certain number of linearly independent, commuting conser- 
vation laws, which are represented by selfadjoint operators Q A , . . . , Qj^ on H A 
commuting with the dynamics of the reservoir, i.e., 

[H A ,Qf] = 0, [Qf,Qf] = 0, [Q A ,a] = 0, (2.2) 

for alH, j, = 1, . . . , M, and for all "observables" a G A A . More precisely, one 
assumes that all operators expitH A /h, {expi Sj Qf} _ x commute with one an- 
other, for arbitrary real values of t,sx,..., s M - 

A typical example of a conservation law is the particle number operator, N A , 
of a reservoir consisting of a gas of non-relativistic atoms. 

On the algebra, B(H A ), of all bounded operators on the Hilbert space H A , we 
define "gauge transformations of the first kind" by setting 

tffa) := e is - QA ae- is - Q \ (2.3) 

for a G B(H A ), where 

M 

s Q A := J2 s i®f- < 2 - 4 ) 

3=1 

Then {ip A | s G M M } is an A/-parameter abelian group of *automorphisms (see 
(I2l6b . below) of the algebra B(H A ), and 

«f (rf(a)) = <p A (at (a)) , (2.5) 

for all a G B(H A ), by (12.21) . It is natural to define the "observable algebra" A A as 
the algebra of all those operators a G B(H A ) for which 

^(a) = a, for all s G R M . (2.6) 

To every conservation law Q A there corresponds a conjugate thermodynamic 
parameter, /i, , commonly called a chemical potential. 

Thermal equilibrium of the reservoir at inverse temperature (5 and chemical 
potentials = (/ii, . . . , /jm) is described by a mixed state, or density matrix, 
given by 

(H^J-'exp-^^-M-QA], (2.7) 

where 

S A jM := tr (exp -(5 [H A - » ■ Q A ] ) (2.8) 
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is the grand-canonical partition function. Of course, it is assumed that the opera- 
tors 

exp -/3 [H A -ii- Q A ] 

are trace-class, for all chemical potentials /x in a region M. C R M , for all (5 > 0, 
and for arbitrary compact subsets A of physical space R 3 . It is also commonly 
assumed that reservoirs are thermodynamically stable, in the sense that the ther- 
modynamic potential, G, given by 

PG&hV) := -ln(H^) (2.9) 

is extensive, i.e., proportional to the volume, V, of the set A, up to boundary 
corrections, for arbitrary (3 > 0, fi G M. 

The expectation value of an operator a G B(H A ) in the equilibrium state 
corresponding to the density matrix (12.71) is given by 

<(«) ■= (^^(e-^-^a). (2.10) 

The state o/g has some remarkable properties to be discussed next. 
It is time-translation invariant, i.e., 

«#mK(«)) = 0)> (2- 11 ) 

for arbitrary a G B(H A ). It obeys the celebrated KMS condition 

<4* («?(«) & ) = <m (^»)) , (2-12) 

for arbitrary a and 6 in B(H A ). In particular, if a G *A A then 

<4# & ) = <m («)) . ( 2 - 13 ) 

for arbitrary b G B(H A ); see (ESI). Eq. (I2T21 is an easy consequence of eqs. (12.101) 
and (12.31) and of the cyclicity of the trace, i.e., 

tr(a5) = tr(6a). 

For further details concerning these standard facts of quantum statistical mechan- 
ics we refer the reader to [Ru 3, BR]. 

Next, we recall some conventional wisdom concerning the thermodynamic 
limit of a reservoir. We are interested in understanding asymptotics of physical 
quantities, as the region A to which the reservoir is confined increases to all of R 3 , 
or to an infinite half- space 

R 3 : = {(x,y,z) G R 3 | x ^ 0} . (2.14) 
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We use the notation "A / oo" to mean that A / M 3 (or A / M|), in the sense 
of Fisher (meaning, in essence, that the ratio between the surface and the volume 
goes to zero); see [Ru 3]. 

We introduce an operator algebra T r , called the "field algebra", convenient 
for the description of the thermodynamic limit of a reservoir: 

F := \j B(H A ), (2.15) 

A /oo 

where Va/oo B(H a ) is the algebra generated by all the operators in the increasing 
sequence of algebras 

... C B(H A ) C B(H A ') C 

A C A', and (■) denotes the closure in the operator norm. Technically speaking, 
T r is a C*-algebra, [BR]. The superscript r stands for "reservoir". Below, we 
consider two interacting reservoirs, labelled by r = /, II. 

A group {r t \ t e K"} of homomorphisms, of a C*-algebra T is an n- 
parameter * automorphism group of T iff 

r t=0 (a) = a, r t (r t '(a)) = T t+ ?(a), 

and 

r t (ay = T t (a*), (2.16) 

for all a 6 T and arbitrary t, t' e R n . 

For the purposes of this paper we shall require the following two assumptions 
concerning the existence of the thermodynamic limit. 

(Al) Existence of the thermodynamic limit of the dynamics and the gauge 
transformations . 

For every operator a G T r , the limits in operator norm 

n — lim af (a) =: a t (a) (2.17) 

A/"oo 

and 

n-lim ip A (a) =: <p a (a) (2.18) 

A /"oo 

exist, for oil f £ 1, s G M. M , and define * automorphism groups of the field alge- 
bra T r . The convergence in \2.17\ and \2.18l is assumed to be uniform for t in 
any compact interval ofM. and for s in any compact subset o/]R M , respectively. 
The * automorphism groups a t and ip s may be assumed to be norm-continuous in 
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t and s, respectively; (but "weak* continuity" will usually be sufficient). 

We define the "kinematical algebra ", A r , to be the largest subalgebra of T r point- 
wise invariant under {y? s }seR M ' i- e -> 

A r := {a G F | <p e (a) = a, for all s G R M } . (2.19) 

Since, by (12.51) . (12.171) and (12.181) . a t and ip s commute, a t (a) G A r , for every 
a G A r . 

(A2) Existence of the thermodynamic limit of the equilibrium state. 

For every a G T r , 

A lim = : wpA a ) ( 2 - 2 °) 

exists and is time-translation invariant, i.e., 

wp t p (a t (a)) = uJp ifl {a), (2.21) 

for all a G J- 7 ", t G R. 

o 

We assume that T r contains a norm-dense subalgebra T r with the property that 
the operator at(<p a (a)) extends to an entire function of t G C, s G C M , for every 

o 

a G T r . If a t and r s are norm-continuous in t and s, respectively, the existence 

o 

of an algebra T r C T T with these properties is an easy theorem. 

From the KMS condition d2Tl2b . and from dTTTl) . (T2~T8l) . (l2!20b . it follows that 
the infinite-volume state u}p )tt on T r obeys the KMS condition 

u)/3 tlx (a t (a)b) = w/j^^fti+p (2.22) 

for all aeT', b G T r . If a G „4 r n / r then eq. (12321 simplifies to 

= W0, M (&a!t+i/Mi(a)). (2.23) 

2.2 Thermal contacts and tunnelling junctions between macro- 
scopic reservoirs 

We consider two reservoirs, I and 77, with all the properties described in Sect. 2.1. 
These reservoirs may or may not have the same physical properties. For example, 
they may be ordinary metals located in two complementary half- spaces of R 3 ; or / 
may be a metal and II a superconductor, etc. Later, we shall consider the example 
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where I and 77 are ordinary metals, i.e., non-interacting electron liquids. In the 
following, "J" will be a shorthand notation for (I, A 1 ), and "77" for (II, A n ), 
where A 1 and A 11 are arbitrary compact subsets of R 3 . Realistically A 7 and A 11 
should not intersect; but we shall ignore this constraint. 

The Hilbert space of the system obtained by composing the two reservoirs is 
given by 

H = H 1 ®H n (2.24) 

and the dynamics, before the reservoirs are brought into contact, is generated by 
the Hamiltonian 

H° := H 1 ®! + 1®H H . (2.25) 

The natural algebra of operators of the coupled system is given by BiTi 1 ) ® 
B(H n ) and, in the thermodynamic limit (A 1 / R 3 _ and A 11 / R\, or A 7 / R 3 

and A 77 / R 3 ), by 

F := T 1 ®? 11 , (2.26) 

where T 1 and T 11 are the field algebras of the two reservoirs (see d2.15l) ). and the 
closure is taken in the operator norm. 

A contact or tunnelling junction between the two reservoirs is described in 
terms of a perturbed Hamiltonian, H, of the coupled system. The operator H has 
the form 

H = H° + W(A I ,A n ) (2.27) 

where W (A 7 , A 77 ) is a bounded, selfadjoint operator on H for each choice of A 7 
and A 77 . We shall always require the following assumption. 

(A3) Existence of the thermodynamic limit of the contact interaction. 

n-lim W (A 7 , A 77 ) =: W (2.28) 

A'/oo 
A 11 foe 

exists (as a selfadjoint operator in T). 

Let af 1 and 1 be the time evolutions of the reservoirs before they are 
brought into contact; see eq. d2.ll) . The time evolution of operators in BiTi 1 ) <g> 
B(7i n ) is then given by ® af u and is generated by the Hamiltonian H° in- 
troduced in (I2.25I) . After the interaction ^(A 7 , A 77 ) has been turned on the time 
evolution of operators in the Heisenberg picture is given by 

a 7u77 ( a ) = e^ tH ^ae- l{tH ' h \ (2.29) 



with H as in (l2~27l for a e B(H T ) <g> B(H n ). 
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It follows from assumptions (Al) and (A3), Eqs. (12.171) and (12.281) . that the 
thermodynamic limit of the time evolution of the coupled reservoirs exists: For 
arbitrary a E T the limits 

a°(a) = n - lim af ® af"(a) (2.30) 

A'/oo 
A'Voo 

and 

a t ( a ) = n - lim a t /u// (a) (2.31) 
k 1 /oo 

A 11 /oo 

exist, and the convergence is uniform on arbitrary compact intervals of the time 
axis. 

Given (12.30b . (12.311) follows by using the Lie-Schwinger series for a t (a°_ t (a)). 
We distinguish between different types of contacts or junctions between the two 
reservoirs, according to symmetry properties of the contact interactions W(A J , A 11 ) 

(Jl) Thermal contacts. The interaction W(A ] , A 11 ) commutes with all the 
conservation laws {Qf ® l}j£i and {11 ® Qf n }j=i of the reservoirs, i.e., 

rf'(W(A',A*)) = <A" (W(A^A")) = W(A',A*), (2.32) 

for arbitrary A 1 and A /7 , where (fg 1 is a shorthand notation for ® id, and 
<^ 7/ stands for id <g> y^"; (see eq. d23b ). It follows from (12.271) and d2~2l that 
the operators {Qj 7 ® 1} and {11 (g> Qf 11 } are conservation laws of the perturbed 
dynamics. 

By assumption (Al), the limits 

ra -lim yf{a) =: y/» (2.33) 

A r /oo 

exist for r = / or 71 and for all a E T . By assumption (A3) and (12.321) it follows 
that 

¥>£W = tf{W) = W } (2.34) 

in the thermodynamic limit. 

Energy appears to be the only thermodynamic quantity that can be exchanged 
through a thermal contact. 

(J2) Tunnelling junctions. There are m < min (M 1 , M 11 ) linear combina- 
tions, Qi' , . . . , and Qf" , . . . , Q^ 1 , of conservation laws of the two reservoirs 
with the property that the operators 

Qf n := Qf ® 1+ 1 ® Qf', (2.35) 
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j = 1, . . . , to, are conservation laws of the perturbed dynamics generated by the 
Hamiltonian H of eq. (I2.27I) . Without loss of generality, we may assume that 
Qf 1 = Qf and (J^ = Qf , for j = 1, . . . , m. Of course, there may be further 
conservation laws of the reservoirs, Qf 1 ®i, i^Qf" , for some i > m and/or some 
j > m, which are conservation laws of the perturbed dynamics. "Leaky junctions" 
are contacts where the interaction W(A I , A 77 ) violates some of the conservation 
laws Qf ® 1 and/or 1 ® Qf , i,j > m. For convenience, we shall sometimes 
assume that the operators Qj 11 , j = 1, . . . , m, are the only conservation laws of 
the perturbed dynamics, and M 1 = M n = m. Let s = (s x , . . . , s m , 0, . . . , 0). We 
define 

^ s UlI (a) := rf'®^), (2.36) 
for a G BiT-L 1 )® BiU 11 ), Mid 

ip a (a) : = n-lim </? 7u77 (a), (2.37) 

A'/oo 
A" /oc 

for a G T\ see assumptions (Al), eq. (12. 1 8b . 

Tunnelling junctions can then be characterized by the requirement that 

tp 1 ^ 11 (W (A 1 , A 11 )) = ^(A 7 ,A 77 ), (2.38) 

for arbitrary A 7 , A 77 , and hence, using (12.181) and (12.281) . we find that, in the ther- 
modynamic limit, 

<Ps(W) = W. (2.39) 

As an initial state of a tunnelling junction we shall usually choose a state uj close to 
a tensor product state, ujfi j £g> uf n n , of two equilibrium states of the uncoupled 
reservoirs, where /3 7 , fj, 1 and /3 77 , /x 77 are arbitrary, (with fi 1 G M 1 , n 11 G M. 11 ). 

Two reservoirs joined by a tunnelling junction can exchange energy and 
"charge" (as measured by the conservation laws Qf 1 ® 1, \®Qf n ,j = 1, . . . , m), 
or leak some "charge" corresponding to ® 1, or to 1 Cg) Qf" , for some j > m. 

Energy current. The operator corresponding to a measurement of the gain 
of internal energy per second of reservoir r, with r = / or 77, at time t is conve- 
niently defined in the Heisenberg picture by 

p r (t) ■.= j t ^ n m 

' a 7u77 ([if,iT]), (2.40) 



h 

where a 7u77 is as in d2~29l . and H r = H 1 ® 11 or = 11 ® if 77 , for r = / or 77, 
respectively. By (12.251) and (12.271) . 

P r (t) = ^a 7u77 ([^(A 7 ,A 77 ),/f r ]). (2.41) 



FMU February 7, 2008 



13 



By (I2.17I) . (12.281) and (12.3 lb . the operator corresponding to the energy gain per 
second of reservoir r has a thermodynamic limit given by 

P r (t) = ~«.KW)U (2-42) 

where a r s is the time evolution of reservoir r in the thermodynamic limit, in the 
absence of any contacts. It follows from (12.411) and (12.421) that 

P\t) + P n {t) = % - a 7u7/ ( [^(A 7 , A 77 ), H°] ) 
= La™([W(A I ,A' I ),H]) 
= -|^(W(A 7 ,A 77 )), (2.43) 
where H° and H are as in eqs. (12.251) . (12.271) . and, in the thermodynamic limit, 

P\t) + P n (t) = -± at (W), (2.44) 

with W e T. 

We observe that if uj is an arbitrary time-translation invariant state of the 
coupled system, we have that 

uo (P 7 (t) + P u {t)) = 0, (2.45) 

for all times. 

Charge current. The operator corresponding to a measurement of the gain 
of charge Q r j per second at time t, in reservoir r, is conveniently defined by 

m = ^<*rm 

= ~ a z t uU ([H,Q^]), (2.46) 



for j = 1, . . . , m, r = I, II, with Q) := Qf <g> 11, and Q 77 := 1 <g> Qf*. Since # 
is given by 

# = H 1 ® 11 + # 77 + 1^(A 7 ,A 77 ), 
and since H 1 <g> 1 and 1 ® H n commute with Q\, it follows that 

m = -^a 7u77 ([^(A 7 ,A 77 ),Q5]). (2.47) 
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In the thermodynamic limit, 

with ip r s as in (12.331) : r = 1,11, j = 1, . . . , M r . Recall that, for j = 1, . . . , m, 
the operators Qj Ul1 = Qf' <E> H + H <S> Qj" are conservation laws of the perturbed 
dynamics, see (12.38b . and therefore 

I-(t)+lf(t) = io£ in ([W(A I ,A n ),Q} m \)=0, (2.49) 

for j = 1, . . . , m; this can be transferred to the thermodynamic limit. Apparently, 
charge lost by one reservoir is gained by the other one. 

For j > m, (12.491 ) does not hold in general and the operators lUt),r = I, II, 
describe the leakage of charge Q r j at the junction. 

2.3 Connections with thermodynamics 

We start by recalling the 1 st and 2 nd law of thermodynamics. 

For the reservoir r, the first and second law of thermodynamics can be sum- 
marized in the equation 

dU Ar = T r dS Ar + /x r • dq Ar -p r dV r , (2.50) 

where U Ar is the expectation value of the Hamiltonian H r in a state of reservoir r 
close to, or in thermal equilibrium, i.e., U Ar is the internal energy of the reservoir 
r; T T the temperature; S Ar the entropy; q Ar is the expectation value of the charge 
Qp j = 1, . . . , M r , in the state describing the reservoir; p r is the pressure, and 
V r = vol(A r ) the volume. The differential "gT indicates that we consider the 
variation of U Ar , S Ar , etc. under small, reversible changes of the state of reservoir 
r (which may include small changes of the region A r ). 

We shall be interested in studying small, slow changes in time of the state of 
reservoirs I and 77, at approximately fixed values of the thermodynamic param- 
eters T r , [f and p r , brought about by opening a contact or junction between the 
two reservoirs. Then U Ar , S A ' , . . . are time-dependent, and (12.501) becomes 

jjAT = T rgAr + . ^ _ p r yr ^ y ^ 

where the "dots" indicate time derivatives. By (12.401) . the energy gain per second, 
U Ar , of the reservoir r is given by 

U^ r (t) = lu IuI1 (P r (t)), (2.52) 
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and the gain in the j charge per second by 

qf(t) = (JJ(t)) , (2.53) 

see (12.461) . where lu IuI1 is the state of the system consisting of the two reservoirs. 
By (12.511) . the change in entropy per second of reservoir r is given by ((3 r := 

1/T) 

S Ar = f3 r (u Ar - ii r ■ q Ar + p r V r ^j . (2.54) 
We define the entropy production rate, £ Iun , by 

£ IUII ;= gA' + £A"_ (255) 

The main property of £ Iun is its sign: thermodynamic systems should exhibit 
positive entropy production, 

£ IuU > 0, (2.56) 

in the limit where A 7 f oo and A 11 /* oo. In [Ru2], Ruelle has proven (12.561) for 
the special case of thermal contacts between infinitely large reservoirs. Below, we 
shall derive (12.561) . under more general conditions, from the positivity of "relative 
entropy"; see also [JP1]. 

We consider a situation in which the state of the system consisting of reservoirs 
I and 77, before a contact or junction is opened, is given by the tensor product of 
two equilibrium states 

u IUII (a) = ^;«wjy(fl), (2.57) 

for any operator a E B(H A ') <S> B(H An ), where the equilibrium states u)j£ „ r 
have been defined in (I2.10I) . The state cu Iun is invariant under the unperturbed 
time evolutions a A r , r = I, II, of the reservoirs. 

At some time t , the contact between the reservoirs is opened, and we are inter- 
ested in the evolution of the state cu Iun under the perturbed time evolution, «[ u7/ , 
introduced in (12.27b . (12. 2% . In particular, we are interested in calculating the rate 
of energy gain, or loss, U Ar [t], the gain or loss of charge j, qf (t),j = 1, . . . , M r , 
per second and the entropy production rate £ IUII (t), under the perturbed time evo- 
lution, in the state uj IuI1 . By eqs. (l2~40l) and (l2~4TT) . 

(j Ar (t) = LU IUU (P r (t)) 

= l^ u// (a[ u7/ ([^(A 7 ,A 7/ ),^ r ])) (2.58) 

and, by (l2~46l 

qf(t) = ^(jT(t)) 

= iu^^tiWiA'^q])). (2.59) 
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By assumption (A2), see (12.20b . the states tu Iun have a thermodynamic limit 

u°(a) = lim oo IUII (a) } (2.60) 
a 1 /oo 

A 11 /oo 

for a E T = T 1 <8> ^ r7/ - It follows from this property, from assumption (A3), and 
from equations (12.3 lb . (12.42b . and (12.48b . that the quantities U Kr (t) and g^ r (^) 
have thermodynamic limits 

V r (t) := lim £7 Ar (t) 

A 1 /oo 
A" /oo 

= -±u°(a t (a:(W)))\ s=Q (2.61) 



and 



jUt) := lim qf{t) 

A 1 /oo ■' 



A 11 /oo 

These limits are uniform on compact intervals of the time axis. By (12.43b and 

7>'(t)+7^(t) = -Iw^otCWO) (2.63) 

and 

Jj{t)+J"(t) = 0, (2.64) 

for j = 1, ... , m. 

Next, we study the entropy production rate for finite reservoirs. Let 

p r := (3^» _1 exp - /F [JT - M r • Q Ar ] (2.65) 

be the density matrix corresponding to the equilibrium state, uijf r , for the reser- 
voir r; see (12.101) . Then 

- In p r = p r [H r - n r ■ Q Ar ] - f3 r G {(3 r , /x r , V r ) ■ 1 , (2.66) 

where the thermodynamic potential G is as in (12.9b . If the confinement region A r 
is kept constant in time, so that V r = 0, then it follows from (12.51b . (12.52b . (12.53b 
and (l2~66ll that 

S Ar = (3 r (lI Ar - ii r ■ if 

= -|u/ u '>r(m/)) 

= ~ tr (p? Vp" a?* QnfT)) . (2.67) 
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Thus, by (l2~55l 



S IUII (t) := 5 IUII (t) 

= — trtfvp *™ (lnpW)) . 



(2.68) 



By (12.611) and (12.621) this quantity has a thermodynamic limit 



m = E /3 r tn*) - M r • ^)] • 



(2.69) 



r=/,ZT 



Integrating (12.681) in time, we find that 



S IUII (t)-S IUII (0)= -tr (p 1 ®p n [a 1 ^ 11 (In p 7 ® p 7/ ) -In p 1 ® p 11 ]) 



This equation shows that S IUII (t) — S IUII (0) is nothing but the relative entropy 
of the density matrix a 7u77 (p 7 (g> p 7/ ) with respect to the density matrix p 1 p u \ 
see e.g. [BR, vol II] for a definition of relative entropy, which differs from ours by 
the sign, and [JP1] for similar, independent considerations. If A is a non-negative 
matrix and B is a strictly positive matrix then 



see Lemma 6.2.21 of [BR, vol II]. Setting A = p 1 ® p" and B = ct /u/J '(p 1 ® p ff ), 
we find that 



S^( t ) _ 5 /u// (0) > ^(p^p^-a^V®/ 1 )) = 0, (2.72) 



by the unitarity of time evolution and the cyclicity of the trace. It follows that 



and this inequality remains obviously valid in the thermodynamic limit: 



(2.70) 



tr (AlnB — A In A) > tr (A - B) 



(2.71) 




(2.73) 




(2.74) 



Thus, if the limit 



lim £{t) =: £ 



(2.75) 



exists then 




(2.76) 
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i.e., the entropy production rate S, in the thermodynamic limit, is non-negative, 
as time t tends to oo; see [Rul]. In Sect. 5, we shall study examples where £ is 
strictly positive. 

Let us assume that the entropy production rate £{t) converges as t — > oo. It 
follows from (12.741) and (12.761) that it is nonnegative. 
Let us set 

V := V 1 = -V 11 (2.77) 

and, for j — 1, . . . ,m, 

Jj ■= J- = -J"- (2.78) 

In the case where each reservoir has precisely m conservation laws {Q^JL^ r = 
J, II (that is, M 1 = M u = m and {Qj <g> t + 1 <g> Qf}f =1 are conservation laws 
of the coupled system), nonnegativity of the entropy production rate implies that 

g = (p 1 - (3 n ) V - (/3 V - /3 7 V 7 ) ■ J > • (2.79) 

The currents J- vanish for thermal contacts, and ( 12.791) shows that energy is trans- 
ferred from the hotter to the colder reservoir — as expected. 

The thermoelectric situation corresponds to M 1 = M 11 = m = 1, and 
Q\ = N r is the particle number operator. For identical temperatures but different 
chemical potentials, ( 12.791) shows that particles are transferred from the reservoir 
with the higher chemical potential to the reservoir with the lower chemical poten- 
tial. Notice that energy may flow from the colder reservoir to the hotter one when 
the chemical potentials are different (consider e.g. (3 1 pJ 3> P n fi n but (3 1 > (3 n ). 

Also interesting is the case of adiabatic thermal contacts between two reser- 
voirs, i.e., without heat exchange. A general discussion of systems with time- 
dependent interactions confined to time-dependent regions is given in [FMSU], 
and will be elaborated upon in a forthcoming paper. 

2.4 Existence of stationary states in the thermodynamic limit 

The above considerations, and in particular (12.611) . (12.621) . and (I2.75I) . suggest to 
study the question whether the infinite- volume states 

uj t (a) := uj (a t (a)), aeF, (2.80) 

have a limit, as t — » oo. 

The state cu°, defined in (I2.60I) . is obviously invariant under the unperturbed 
time evolution defined in (12. 301) . Thus 



(2.81) 
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A sufficient condition for the existence of a stationary (i.e., time-translation in- 
variant) limiting state, 

^stat(a) = lim u t (a), a G T , (2.82) 

t— >oo 

is given in 

(A4) Existence of a scattering endomorphism. 

The limits 

o±(a) = n — lim a°_ t (a t (a)) (2.83) 

t— +dzoo 

exist, for all a G JF, and define * endomorphisms of J 7 , i.e., a± are homomorphisms 
of the C* -algebra T with the property that cr ± (a)* = cr±(a*),for all a G T. 

The usefulness of these so-called scattering (or M0ller) endomorphisms has first 
been recognized in [He, Rob]; interesting examples have been constructed in 
[BM]. In the context of thermal contacts and tunnelling junctions, they have first 
been used in [DFG]; see also [Ma]. 

It is important to note that scattering endomorphisms do not exist in finite 
volume, because the free and the perturbed time evolutions of the two reservoirs 
are generated by Hamiltonians 

H° = H 1 <g> t + 1 <g> H", 
H = H° + W{A I ,A n ), 

see (12.251) and (12.271) . with pure-point spectra when A 1 and A 11 are compact. It is 
thus natural to wonder about the meaning of scattering endomorphisms for large 
but finite reservoirs. Let us sketch the answer to this question. We fix an arbitrarily 
small, but positive number e. For every operator a G T, there exist compact 
regions A r (e, a), r = 1,11, and an operator a £ G B(H A ) (g> B(H A ), with A r = 
A r (e, a),r — I, II, such that 

e 

\\a — a £ \\ < — . 
n en 4 

Then, by (12.821) and (12.831) . there is some T{e, a) < oo such that 

l^statO) - w t (a £ )\ < I , 

for all t > T(e, a). Assumption (A3), eq. (12.311) . tells us that, for an arbitrary T < 
oo, there are compact sets A r (e, a, T) D A r (e, a) such that if A r D A r (e, a, T), 
r = I, II, then 

|K(a £ ) -a[ u// (a £ )|| < | , 
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for all t G [0, T] . Finally, by assumption (A2), one can choose A r (e, a, T) so large 
that 

|^K^(a £ ))-^°W UJ/ (0)| < f , 

provided A r D A r (e, a,T), r — I, II, for all times t G [0, T). It follows that, for 

any T, with < T{e,a) < T < oo,andforA r D A r (e,a,T), r = I, II, 

\u stat (a) - u Iun (a[ u// (a £ )) | < e , (2.84) 

for all times t, with T(e,a) <t < T. 

These simple considerations, combined with (12.441) and d2.48l> . show that the 
energy-gain rates U Ar (t) and the currents q 1 ^ (t) of two very large, but finite reser- 
voirs, r = I, II, are well approximated by the energy-gain rates 

V r := Km V r (t) = -±u; stat (a:(W))\ s=0 (2.85) 

and the currents 

JJ := lim JJ{t) = ~~ ^tat mW)) | s ^ , (2.86) 

respectively, for a large range of sufficiently large, but not exceedingly large times 
t; (see (l2~6T1) . (12321 ). 

Remark. It is usually much easier to prove that the limits 

&'±(a) := n — lim ct- t (a°(a)) (2.87) 

exist and are operators in T , for arbitrary a E T , rather than to establish the exis- 
tence of the scattering endomorphisms a± in (12.83ft . If the unperturbed dynamics 
of the reservoirs is dispersive, (as for non-interacting, non-relativistic electrons), 
one may hope to prove (12.871) by using a simple Cook argument; see e.g. [He, 
Rob, CFKS]. If both limits (l2~83b and (l2~87l) exist then 

<r±K(a)) = <4M<0) = ( 2 - 88 ) 

i.e., cr^_ is a left and right inverse of a±, and hence cr± is a * automorphism of 
This will turn out to hold in the examples discussed in subsequent sections. 

2.5 Uniqueness and stability properties of stationary states 

We first describe the property of return to equilibrium for a single reservoir. Let 
u be a state on the field algebra of a single reservoir, T r , i.e., u is a positive, 
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linear functional on jF r normalized such that uj(1) = 1. From jF r and uj one can 
construct a Hilbert space H w , a representation tt^ of jF r on H^, and a unit vector 
Vt E ri w (unique up to a phase) such that 



= {n w (a)n | a6?}, (2.89) 

where the closure is taken in the norm on ri u , i.e., is "cyclic" for n u} (J rr ), and 

u(a) = (fi,7r w (a)fi), (2.90) 

where (•, •) is the scalar product on Ti^. This is the content of the Gel'fand- 
Naimark-Segal (GNS) construction. If the state uj is time-translation invariant 
then there exists a one-parameter unitary group {U u (t) I i 6 R} on ri w such that 

7r w (04(a)) = C4,(0Tw(fl)Cw(*)*, 

and 

t/ u (t)fi = fi. (2.91) 

Under standard continuity assumptions on U w {€), we can summon Stone's theo- 
rem to conclude that 

U u (t) = e itL " /h , (2.92) 

where the generator L u is a selfadjoint operator on H u with L u Vt = 0. 

A state p on T r is called normal relative to uj iff there exists a density matrix, 
P, on such that 

p(a) = ti Hu (P7r w (a)), (2.93) 

for all a e JF*. 

Let c<j := Wjg )M be an infinite- volume equilibrium state on T r obeying the 
KMS condition (12.22b . We assume that the Hilbert space obtained from the 
GNS construction is separable and that the cyclic vector Q, G H w is the only 
eigenvector (up to phases) of the operator L w of ( 12.921) . which, in this context, is 
called the Liouvillian or thermal Hamiltonian. In other words, the spectrum of L u 
is purely continuous, except for a simple eigenvalue at 0. 

This assumption implies the property of "return to equilibrium" : If p is an 
arbitrary state normal relative to uj = w ft)X then 



lim - / dtp (ot(a)) = u(a), (2.94) 

T-*oo 1 J Q 

for all a e T r . If the spectrum of L w is absolutely continuous, except for a simple 
eigenvalue at 0, then 

lim p (04(a)) = Lu(a), (2.95) 

t— >oo 
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for all a E T r . Eqs. (12.941 ) and (12.951) follow from our assumptions on the spec- 
trum of L u and the KMS condition (12.221) : see e.g. [BFS]. 

Assuming the existence of the endomorphism a + , see (A4), we now address 
the question of uniqueness and dynamical stability of the stationary state a; sta t m 
(l2~82l 

We suppose that the property of return to equilibrium, ( 12.951 ). holds for each 
reservoir separately; (this can be shown for reservoirs consisting of free fermions 
as considered in Sections 3 and 4). Recalling that the reference state u; is a 
product of two KMS states (see (12.601) . (I2.57I) '). it is not difficult to extend the ar- 
guments in Sect. Ill, D of [BFS] to show that if p is an arbitrary state on T normal 
relative to the state tu° in (12.601) then 

lim p(a° t (a)) = ui°(a), a E T , (2.96) 

where a° is the time evolution of the reservoirs before they are coupled; see (I2.30I) . 

Eq. ( 12.961 ) and the existence of a scattering endomorphism cr + , see eq. (12.831) . 
now imply that 

lim pt{a) = lim p(a t (a)) 

t— >oo t— >oo 

= lim p(a° t (a + {a))) 

t— >oo 

= w°Ma)) 

= ^stat(a) , (2.97) 

for a E T . Thus, if the initial state p is an arbitrary state normal relative to cu° 
then the states p t tend to the stationary state a; sta t, as t — > oo, (uniqueness). 

Next, let p be an arbitrary state normal relative to cj stat . We claim that if 
eq. (12.951) holds for each reservoir then 

lim p t (a) = ^stat(a), a E T. (2.98) 

Eq. (12.981) is the property of "return to the stationary state" (stability), for states p 
normal relative to o; s t a t . To prove eq. (12.981) . we follow the arguments in Sect. Ill, D 
of [BFS]: Since p is normal relative to u; s tat, there exist non-negative numbers, 
p n , n = 1,2,3, .. ., with Y^=iPn = 1. and nets of operators {<} aG /„, n = 
1, 2, 3, . . ., with E T, for all a and all n, such that 

oo 

P(«) =y2Pn lim <7r(<)fi,7r(a)vr(<')fi) 

' * a, a' 

n=l 

where n is the GNS representation and VL the cyclic vector corresponding to 
(^stat, F)\ see d2~89l) . (12^01) . and (HH). Then 

oo 

Pt(a) = y^Pn lim <vr(<)fi, 7r(a t (a))7r(<')fi) . 

* ql.ol' 
n=l 
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Let 

6i := < , & 2 := <' , 
for some fixed n, a, a'. Then 

(tt^Q, 7f(a t (a))n(b 2 )Q) = uj stat (b{a t (a)b 2 ) . (2.99) 

Since uj stait (a) = uj (a + (a)) . and since, by eg. (12.831) . 

<7+(a t (a)) = a?(o+(a)) , (2.100) 

for arbitrary a & J 7 , the R.S. of (12.991) is given by 

Wst»t (6*04(0)62) = w°(o- + (6*) a£(o+(a)) o" + (6 2 )) . 

It follows from (12.961) by polarization that 

Km u°(a + (bt)a° t (a + (a))a + (b 2 )) = uj\a + (b\) a+(b 2 )) oo°(a + (a)) 

= u; sta t(6* 6 2 ) ^stat(o) 

= <7r(6i)fi,7r(6 2 )fi) ^stat(a) . (2.101) 

Our contention, eq. (I2.98I) . follows from this. 

2.6 Cluster properties and profiles in d>2 dimensional 

systems 

The last question we wish to address, in this summary of the general theory, con- 
cerns cluster properties of the stationary state u; s tat, which will show that co> stat 
cannot be an equilibrium (KMS) state for the dynamics, a t , of the coupled reser- 
voirs and is, in general, not normal relative to the product state, cu°, of the uncou- 
pled reservoirs. 

We consider two increasing families of reservoirs confined to regions A r C 
R 3 , with 

A r / R 3 , r = IJI, (2.102) 

joined together by a thermal contact or a tunnelling junction localized near the 
origin, x = 0, of physical space. The more realistic situation where the reservoirs 
are confined to two complementary half spaces, R 3 and R 3 , respectively, with a 
junction localized near the origin, has been considered in [DFG]; see also [Rul, 
Ru2]. It will be studied in more detail elsewhere. In order to describe spatial 
properties of the system, we make the following assumption. 

(A5) Existence of space translations. 
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For each reservoir r — I, II, there exists a * automorphism (semi-) group 

{t£ |*GM? ±) } (2.103) 

of the field algebra T r , representing space translations o/IR 3 (Wt±, respectively) 
on T r . 

For the system of two coupled reservoirs, 

r x := t x ®t x ! , xeR 3 , (2.104) 

defines a representation of space translations as a 3-parameter group of *auto- 
morphisms on the field algebra T = T 1 ® J 711 . It is plausible that space transla- 
tions satisfy the following assumption. 

(A6) Asymptotic abelianness of space translations, and homogeneity 
of reservoirs. 

The action of r x on T is norm-continuous in x 6 M 3 and for all operators a 
and b in T, 

lim || [r x (a),b] \\ = . (2.105) 

Furthermore, the dynamics and the equilibrium states of the uncoupled reservoirs 
are homogeneous, in the sense that 

a° t (r x (a)) = r x (a° t (a)) (2.106) 

and 

Lu°(r x (a)) = Lu°(a) , (2.107) 

for all a G T. 

The local nature of the perturbation, W, of the dynamics of the system due to the 
contact or junction, see Assumption (A3), eqs. (12.271) and (12.28b . and assumption 
(A6), then imply that, for all a E T, 

lim \\a t (r x (a)) - a%r x {a)) \\ = 0, (2.108) 

| cc I — >oo 

for all times t. A proof of (12.1081) follows from the Lie-Schwinger series for 
a_ t («°(a)) and use of (12.281) and (I2.105I) . Relation (12.1081) shows that observables 
localized far from the junction evolve according to the non-interacting dynamics. 

It is tempting, and can be justified in examples, to strengthen assumption (A4) 
(existence of scattering endomorphism) as follows. 



FMU February 7, 2008 



25 



(A7) Cluster properties of the scattering endomorphism. 

The limits 

n-lim a°_ t (a t (r x (a))) = a ± (r x {a)) (2.109) 

t— >±oo 

are uniform in x G M. 3 , for every a G T. 

Equations (12.1081) and (12.1091) imply that 

lim \\a ± (T x (a)) - T x (a)\\ = 0, (2.110) 

for every a G T. From this property we conclude that 

lim uj stat (r x (a)) = lim u° (a + (r x (a))) 

|cc|^oo |a;|^oo 

= lim LU (r x {a)) 

\x\^oo 

= co°(a), aeJ 7 , (2.111) 

i.e., very far from the junction, the stationary state c<j stat resembles the product 
state Ld° of the uncoupled reservoirs. 



Remark. It is not hard to understand that if the two reservoirs occupy com- 
plementary half spaces, and Mi, then (12 .111b is replaced by 

lim w s tat (r (a; ,o,o) (a ® ll)) = uj^i i(a) , 

for a G J 71 , and 

for b G T n . This may prove the presence of a profile of temperature, density, or 
... in the stationary state, u; sta t 5 of the system. 

In the examples studied in subsequent sections, which concern junctions be- 
tween ordinary three-dimensional metals, assumptions (A6) and (A7) can be ver- 
ified. Instead of three-dimensional reservoirs, we could consider one-dimensional 
(wires), or two-dimensional (layers) reservoirs joined by a thermal contact or a 
tunnelling junction. Our analysis in Sect. 4 will show that, in dimensions d = 1,2, 
assumption (A7) will fail, in general. In fact, the example of the one-dimensional 
XY spin chain treated in [DFG, AP] and the example of quantum wires studied 
in [ACF] show that, in one dimension, c<j stat may well be space-translation invari- 
ant, i.e., it does not exhibit any profile. In the example of quantum wires, to sta ,t is 
actually a homogeneous thermal equilibrium state. 
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Thus, we observe that the validity of assumption (A7) critically depends on the 
dimension of the reservoirs. In dimension d > 2, this assumption can be expected 
to hold, while it usually fails in dimension d = 1,2. For people familiar with 
elementary facts of scattering theory this will not come as a surprise. 

In a subsequent paper, we will show that, for a large class of reservoirs, one 
can construct "observables at infinity" , see e.g. [BR, Vol. II], corresponding to 
the operators P r \t) and lUt) defined in eqs. 92.421 . (12.481) . respectively. Clearly, 
the expectation values of these operators vanish in the product state u; of the 
uncoupled reservoirs and are given by V r and J J in the stationary state, cu stat , of 
the coupled reservoirs. If we can show that V r ^ 0, or that JJ ^ 0, for r = I 
or II and some j, then it follows that LO stsA is not normal relative to cu°. In the 
examples studied in Sects. 4 and 5, we shall encounter instances where V and 
J T do not vanish. 



3 Reservoirs of non-interacting fermions 

This section serves to introduce a class of simple, but physically important exam- 
ples of reservoirs to which the general theory outlined in Sect. |2] can and will 
be applied. Our examples describe a quantum liquid of non-interacting, non- 
relativistic electrons in a normal metal or a semi-conductor, possibly subject to an 
external magnetic field, or an ideal quantum gas of fermionic atoms or molecules. 
In a subsequent paper, we shall also consider examples describing chiral Luttinger 
liquids, which arise in connection with the quantum Hall effect. 

We start by considering a system consisting of a single, non-relativistic quant- 
um-mechanical particle confined to a region A of physical space M. d , d = 1, 2, 3. 
The Hilbert space of pure state vectors of this system is given by the space 

L 2 (A,d d x) (3.1) 

of square-integrable wave functions with support in A. If the particle has spin 
and/or if there are several species of such particles then L 2 (A, d d x) must be re- 
placed by the space 

h A := L 2 (A,d d x) <g>C fc , (3.2) 

where k = X^iO^a + 1)> is the spin of species a, and I is the number of 
species. 

The one-particle dynamics is generated by the following selfadjoint operator, 
t A , acting on h A , 
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where M is the mass of the particle and A is the Laplace operator on L 2 (A, d d x) 
with selfadjoint boundary conditions (e.g. Dirichlet, Neumann, or periodic) im- 
posed at the boundary, <9A, of A. 

In the following, we choose units in which h = 1 and M = |. 

Other operators are physically interesting. Electrons in semi-conductors would 
involve a potential operator that is diagonal in the space representation. A mag- 
netic field could also be considered; the Laplacian should be replaced by the co- 
variant Laplacian, and a coupling between the spin of the particle and the magnetic 
field should be introduced. In this paper, we restrict our attention to the situation 
<l3~3l . 

Next, we consider a system consisting of n identical particles of the kind just 
considered, all confined to the region A. Its state space is given by a subspace of 
the n-fold tensor product of h A of fixed symmetry type, 

hi := P(h A )® n , h K := C, (3.4) 

where P is the orthogonal projection onto the subspace of wave functions of the 
selected symmetry type under permutations of the n particle variables. If the 
particles are bosons then P = P + projects onto completely symmetric n-particle 
wave functions; while, for fermions, P = P_ projects onto totally anti- symmetric 
wave functions. In this paper, we focus our attention on fermions. 

If the particles do not interact with each other the Hamiltonian, T A of the 
n-particle system is given by 

n 

T n A := ^11®...®^®...®!, (3.5) 

i=i 

where t A acts on the j th factor in the n-fold tensor product in (13 .4b . 

If the number of particles can fluctuate (e.g., because the system is coupled to 
a particle reservoir such as a battery) then it is convenient to use the formalism of 
"second quantization", which we briefly recall. 

The Fock space is defined by 

oo 

H A := hi . (3.6) 

n=0 

The free dynamics on H A is generated by the Hamiltonian 



oo 
n=0 



(3.7) 
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with T A as in (13.51) . The particle number operator, N A , is defined by 

oo 
n=0 

Let k be a symmetric fc x k matrix acting on C fe . We set 

n 

K A : = 1 ® . . . <g> (l <g> Kj) <g> . . . <g> 1 , 
j'=i 

where acts on the j th factor, ~ L 2 (A, d d x)®C k , in the n-fold tensor product 
defining h A . For t A as in (I3.3I) . a typical charge operator, Q A = Q a (k), is of the 
form 

oo 

Q A ■= . (3.9) 

n=0 

The operators H A ,N A and Q a (k) are unbounded, selfadjoint operators on 
ft A ; see e.g. [RS]. 

Next, we describe the structure of H A in some more detail and introduce cre- 
ation and annihilation operators. Let x,y, . . . denote points in physical space M. d , 
and let s = 1, . . . , k label on orthonormal basis in C k . Vectors /„ in the n-particle 
space h A can be represented as square-integrable wave functions, 

fn (^1) • • • 5 VCni ^n) i 

with support in A dn C R. dn , which, for fermions, are totally anti-symmetric under 
permutations of their n arguments. Vectors ■0, 0, ... in Fock space correspond to 
sequences, 

1> = {f n )Zo , = (<?X=o,--- (3-10) 
of n-particle wave functions in h A . The scalar product on H A is defined by 

oo „ n 

(x h st,...,x n , s n ) g n (xi, si, . . . , x n , s n ) . 

n=0 si,...,s n j=l 

(3.11) 

The vector represented by the sequence (/ n )^=0' w ^ m fo — 1> fn = 0, for n > 1, 
is denoted by and is called the vacuum (vector). 

Let V = V A be the linear domain of vectors ip — (fn)n°=o m ^ A w i m me 
property that all but finitely many f n 's vanish. Clearly, V is dense in H A . For 
/ G h A , we define an annihilation operator, a(f), by 

(o(/)V , ) n (a'i>si J --->a5n,s n ) := 

+ 1 52 / dse f(x,s)f n+1 (x, S,Xi,S 1 ,...,X n ,8 n ) , (3.12) 
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for arbitrary if) = (f n )%Li G and 

o(/)0 := 0. (3.13) 

The creation operator, a*(f), is defined to be the adjoint of a(f) on H A and is 
easily seen to be well defined on V. 

It is well known, see e.g. [RS, BR], that, for fermions, the following "canoni- 
cal anti-commutation relations " (CAR) hold: 

{a*(f), a#(g)} = 0, (3.14) 

for arbitrary /, g in h A where a # = a or a*, and {A, B} := AB + BA is the 
anti-commutator of two operators A and B; 

{a(f), a*(g)} = (f,g) ■ 1 , (3.15) 



where (/, g) := J2 S f\ dx f(x, s) g(x, s) is the scalar product on h A . For bosons, 
(13.141) and (13.151) hold if anti-commutators are replaced by commutators (CCR). 
Formally, 

= ^2 f A dx f( x ' s ) 



and 



with 



l *(f) =^2 dxa*(x,s)f(x,s) 



{a(x, s), a*(x', s')} = 6 SS > 5 {d \x - x') . (3.16) 



A remarkable consequence of the CAR is that the operators a(f) and a*(f) 
are bounded in norm by 



Hf)\\ = h*(f)\\ = 11/11 := y/UJ)- (3-17) 
To see this, we choose an arbitrary ip G V and note that 

||a(/)V|| 2 + \\a*(fW 

= (a(f)ij, a(f)ij) + (a*(/)V>, 
= (^{a(/),a*(/)}V> 
= (/,/)MV>>, 

so that 

< 11/11 ■ HV'II • (3-18) 
Equality in (13.181) is seen from examples. 
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Eq. (13.171) is false for bosons, a(f) and a*(f) being unbounded operators. 

For fermions, polynomials in a(f), a*(f), f € h A , form a *algebra of opera- 
tors on H A which is weakly dense in B(H A ). The "observable algebra" A A is the 
norm closure of the algebra of these polynomials in a(f), a*(f), f E h A , which 
commute with the number operator iV A and, possibly, with further charge oper- 
ators Q a (k), for certain choices of k. Every monomial in a and a* belonging to 
A A has equally many factors of a and a*, since it must conserve the total particle 
number. A general monomial in a and a* is Wick-ordered if all a*'s are to the left 
of all a's. 

In terms of creation and annihilation operators, the operators H A , N A and Q A 
can be expressed as follows. 

H A = V / dxa*(x,s)(t A a)(x,s) , (3.19) 

A^ A = V / dxa*(x,s)a(x,s) , (3.20) 

and 

Q a (k) =2_. / dx a*(x, s) k ss i a(x, s') . (3.21) 

s,s> 

In the examples discussed below and in Sects. 4 and 5, we usually regard 
iV A = Q a (k = 1) to be the only conservation law, besides H A , relevant for the 
description of the reservoirs. In a general discussion, we consider M conservation 
laws, Q A = Q a (kj), j = 1, . . . , M, and choose t A as in (13.3b . 

The main result of this section is the following theorem. 

Theorem 3.1. Fort A asin A3. 31 . andQ A = Q A (Kj),j = 1, . . . , M, with /c 1( . . . , k m 
arbitrary, commuting symmetric kx k matrices, the equilibrium states cu A intro- 
duced in eq. (12.701) exist, for arbitrary j3 > and /J, G M M . 

Assumptions (Al), (A2), (A5) and (A6) of Sect. 2, concerning the existence of 
the thermodynamic limit, A /* M. d , hold. 

The proof of Theorem 13. II is standard. A careful exposition can be found in 
[BR], Sect. 5.2. 

In Section 4 we shall consider a system consisting of two identical reservoirs, 
I and 71, both composed of non-interacting, non-relativistic fermions confined 
to some region A = A 7 = A 11 of M d . A convenient notation for creation and 
annihilation operators for the two reservoirs is the following one. 

a*(x,s,I) := a*(x, s)\ nAl <g> l\ HA a , 
a*(x,s,U) := l| wA i ® a*(x,s)\ HAlI . (3.22) 
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We note that all the operators a # (/, /) = a* (/) ® 1 commute with all the operators 
a*(/, il) = 1 ® If, for convenience, we prefer that they anti-commute we 

can accomplish this feature by a standard Klein-Jordan-Wigner transformation: 

a*(f,I) » a*(f,I), 
a*(f,II) i-> a*(f,n)e^ NAlm K 

The operators on the R.S. of (13.231) will again be denoted by a # (/, r), r 
We introduce the following notation. 

X := (x,s,r) e R d x {1,..., k} x {/,//}, 
^ := (X U ...,X N ) , 



(3.23) 
J, II. 



(JV) 



(si, 



,r^)G{/,/J}^ 



E E 

N 

n 



r=7,/J s=l 
N 



dx. 



dX 



3i 



iV 



(3.24) 
(3.25) 
(3.26) 
(3.27) 
(3.28) 

(3.29) 
(3.30) 
(3.31) 



with a* = a* or a. 

We are now prepared to describe the interactions, jy(A 7 , A 11 ), (see eq. 12.271) ). 
corresponding to thermal contacts or tunnelling junctions between the two reser- 
voirs. We shall always assume that the total particle number of the system consist- 
ing of the two reservoirs is conserved. Thus the interaction Hamiltonian W(A T , A 11 ) 
must commute with the operator 



N iun ._ n a1 (g)l + l(g) N 1 ' 
(see eq. (12.351) ). It follows that W(A r , A 11 ) must have the form 



W(A I , A n ) =J2W N (A I ,A 11 ) 



(3.32) 



(3.33) 



N=l 



where 

W N (A J ,A n ) 



dx (N) / rfyWa*^))^^^^)^), 
Ja n 

(3.34) 



FMU February 7, 2008 



32 



and, for each choice of s 1; r 1; . . ., s^, r N , s\, r[, . . ., s N , r' N , 

( (x™ ,s^ N \r^), (y W , s' {N) , r' (A °) ) 



A? 



is a smooth function of x^ G A N and e vanishing if x^ A N or 
$ A N . In Sect. 4, we introduce weighted Sobolev spaces, Wat equipped with 
norms || | (•) II I jv with the property that 

U^aKAV)!! < IlK'^llk, (3-35) 

for all N = 1, 2, 3, . . .. We shall assume that, for each N, there is a function 
wn £ Wat such that 



9{w) 



/] IH^jvIIIjv < oo j (3.36) 



JV=1 



where u> = (wn)n=i, and 



oo 

U ah 



lim y^|||w^' A — watIIIat = 0. (3.37) 



AV°° 



It then follows that 



exists, and 



ri -lim W^A^A 11 ) =: W (3.38) 

k 1 /oo 
tJ 1 /oo 



\\W\\ < g{w) . (3.39) 

Thus, assumption (A4), eq. (l2~29l . of Sect. 2 follows from (13341 . (l3~35t . (l3~36ll 
and (l3~37l 

If we want to describe thermal contacts we shall require that 

[iV A '®l, iy(A 7 ,A 7/ )] = [1 <g> N A " ,W(A J ,A n )} = 0, (3.40) 

while, for tunnelling junctions, only 

[N IL>11 ^(A 1 ,A n )] = (3.41) 

is required, for arbitrary A 1 = A 11 = Ac M d . 

To conclude this section, we remark that a system of two reservoirs of non- 
interacting fermions, with a one-particle Hamiltonian t A as in Eq. (13.31) . and with 
interactions W{A I ,A n ) as in (I3.34I) - (I3.38I) . satisfies Assumptions (A1)-(A3) and 
(A5), (A6) of Section 2. Assumptions (A4) and (A7) are established in the next 
section for d > 3 and small g(w). 
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4 Existence of M0ller endomorphisms 

The goal of this section is to illustrate the general theory of Section|2]by providing 
a complete mathematical description of a concrete system, namely two coupled 
reservoirs of free fermions in dimension d > 3. An illustration can be found in 
Fig. [H The reservoirs are infinite and without boundary, and the coupling W is 
localized near the origin, in the the sense that ||W||', or ||W||", is finite (see (14.3b . 

d4]§, or 623 ). 



Figure 1: Two-dimensional reservoirs coupled by a local interaction W . We actu- 
ally consider the three-dimensional analogue of this situation. 

The Hamiltonian for each reservoir has been introduced in Section|3l see equa- 
tion (I3.7I) . The coupling between reservoirs is represented by the interaction in 
equations (13.331) and (I3.34I) . As stated in Theorem l3.ll there exist time evolution 
automorphisms a° and a t in the thermodynamic limit. The former corresponds to 
the free dynamics and the latter to the dynamics for interacting reservoirs. 

In this section, we establish the existence of the M0ller endomorphisms a± 
defined in (12.83ft . We start with the Dyson series for a°_ t a t , namely 



where we set W(t) : = a°_ t (W). Convergence of this series for finite t is clear 
since ||W|| < oo, but we need to consider the limit t — > oo. Let us define an 
operator D m (t) on the field algebra T by 





dt x ...dt m [W(t m ),...,[W(t 1 ),a]...], (4.1) 




(4.2) 



for arbitrary a G T . It is understood that D a = a. 
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We define the norm of an interaction W by setting 

k k 

\W\\> = ]T N*W J2 S 

N>1 Sl ,...,s N =l s' v ...,s' N =l 



E £ IK((-,5W,rW),(, S 



/W /(2V) 

1 ' 



n,...,r N =I,II r[, 



the function u>jv is viewed above as a function on M 2dJV , 
defined by 



(4.3) 



and the norm || ■ \\' M is 



23M/2 



dx {M) f(x( M ))H 



k=l 



d 2 



dx ^ 



+ 4 + 1) /(*(">) 



1/2 



• (4.4) 



Note that the operators 



. d 2 
dx'l 



+ x\ + 1 are bounded below by 2, which implies the 



inequality (13.351) . namely ||/||l2 



< 



' M ; this inequality is saturated when / 



is a product of Gaussians centered at the origin, f(x^ M '] 
Theorem 4.1. For d>3we have the bound 

\\[W(t),D m ^(t)a]\\dt<^[-^ 



T-fM 



a||' 



-4/2 



A similar statement holds for negative times. That is, one can rewrite the 
Dyson series in (14.11) and the operator DM(t) in (14.21) for t < by integrating over 
negative times > t\ > • ■ ■ > t m > t. Then Theorem 14. 1 1 holds with an integral 
from — oo to 0. The proof for negative times is identical to the one of positive 
times. 

Before proving Theorem 14. 11 let us work out its main consequence, the exis- 
tence of M0ller endomorphisms. 

Corollary 4.2. If^\\W ||' < 1, there exist o ± such that 

lim \\a°_ t a t (a) — cr±(a) \\ = 

t— >±oo 

for all a with \\a\\' < oo. 

Corollarv l4.2l implies the existence of a scattering automorphism, see Assump- 
tion (A4). We comment below that the norm (14.41) can be replaced by an object 
that is translation invariant, see (14.251) . Therefore the scattering automorphism is 
given by a limit of infinite times, and this limit exists in norm, uniformly with 
respect to space translations (for both the interaction and the operator a E JF). 
Hence Assumption (A7) holds. 
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Proof of Corollarv \4. 21 Observe that 



a°_ t a t {a 



) = a+ Y^i m [ t [W(s),D m _ 1 (s)a]ds. (4.5) 

m>l 



Then, for t < t', 

\\a°_ t a t (a) — a°_ t ,a t i(a) 



< Yl / \\[W(s),D m ^(s)a]\\ds. 

m>l Jt 



By Theorem 14. II and the dominated convergence theorem, the right side vanishes 
as t, t' — > oo. This implies the norm-convergence of a°_ t a t (a). □ 



We will make use of Hermite functions in the proof of Theorem 14.11 so we 
collect a few useful facts on them. The Hermite functions are denoted by {0 9 } 9G n, 
where 

v ^T 7r 4 \dxJ 
with x E R. These functions satisfy the equation 



'2 



(-^2 + ar 2 )^(») = (2g + (4.6) 
Lemma 4.3. 

(nj |(e^0 p ,0 9 )|<^M^. 

a/ 47T|E| 

Proo/ By Cauchy-Schwarz and since the operator — -£3 is positive definite, we 
have that 

1 



y|i= / Vx 2 + 1 \(j) q (x)\ — =dx 

yx 2 + 1 

, d 2 9 \ \i/2/ Z" 00 \V2 



The first factor is equal to >j2(q + 1) and the second one equals y/2n, which 
proves (i). Claim (ii) immediately follows from 



□ 



FMU February 7, 2008 



36 



Hermite functions form an orthonormal basis of L 2 (M). We use them to ex- 
press the interaction as a polynomial of creation and annihilation operators of 
fermions in states described by Hermite functions. The free time evolution of the 
interaction can be described as an evolution of these functions, and their decorre- 
lation in time can be controlled using Lemma l4~3l (ID. Finally, Hermite functions 
will be removed at the expense of introducing differential operators in the defini- 
tion of the norm of the interaction; see (14.41) . 

We use from now on the following notation: for q = (gi, . . . , q^) G N d and 



x = (x u ...,x d ) e R d , 



8=1 



Proof of Theorem W7\ We start by rewriting the interaction in the basis of Hermite 
functions. Let := (Qi, . . . , Qn), with 

Qj = (qp SJ, rj) € n d x {1, . . . , k} x {I, II}. (4.9) 

We set 

r r N 

(4.10) 

Here, X^ N ' is determined by and x^ N \ namely Xj = (xj, Sj, r j 

) if Qj = 

(qr -, Sj, Tj). The interaction (13.321) . (13.331) is given by a sum over 

Q:= (N,QW,Q'W), 

Q Q 

where a* (Q^ ) = Ylf=i a * (Qj ) > an d a * (Qj ) is tne creation operator for a fermion 
in the state (f> q , of spin Sj, in the reservoir rj. The annihilation operators &(Q'^) 
are defined similarly. The operator a also has a Hermite expansion 

a = J2 a Q = £a*(QW)MQ ( ^^W (A °)- (4-12) 

Q Q 

With this notation, we have that 

/*oo P 

/ ||[W(0,^ m -i(t)a]||^< / dt l ...dt m 

JO D D J<X»t m >" >tl>0 

\\[W Q Jt m ),...,[W Ql (t 1 ),a Qo )...)\\. (4.13) 
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The multiple commutator above involves operators W and a, which in turn in- 
volve creation and annihilation operators of particles in both reservoirs. The latter 
satisfy anticommutation relations for particles in the same reservoir, or commuta- 
tion relations for particles in different reservoirs. This introduces a complication 
when estimating the multicommutator above. This complication can be avoided 
by using the Klein- Jordan-Wigner transformation explained in (13.23ft . For sim- 
plicity, we keep the same notation, but we assume from now on that all creation 
and annihilation operators satisfy anticommutation relations. 

Because a°_ t is a *-automorphism, its action on the interaction W simply 
amounts to replacing operators a#(Q) = a#(0 9 , s, r) by 

a°> # (0 q ,s,r)) = a # (e 4 '%,s,r) :=a # (Q,t). (4.14) 

We note that Lemma l4~3l vields the bound 

\\{a _ t (a#(Q)),a _ t ,(a*(Q'))}\\ = |(e^ A q , e*'%OI<W<W 

d 

< (lA^) d/2 nfe + 1 ) l fe' + 1 )^ ( 4 -!5) 
i=i 

The multicommutator in (14.131) can be written as 

[W Q Jt m ),...,[W Ql (t 1 ),a Qo ]...} = 

m 
3=1 

[a*(Qf l) ,t 1 )a(Qi^\t x ),a*(gS Vo) )a(g' w )] ...]. (4-16) 

Here, we set 

a#(^ ) ,t i ) = n«-* J (a # (^)). ( 4 - 17 ) 

where Qjj is the £-th element of Qj. 

A commutator of products of operators can be expanded according to contrac- 
tion schemes. The following equation holds when ki is even: 

[ai . . . a k , h . . . bt] = • • • Oi-ih . . . 

l<i<fc 

i<j<e 

{ai, bj}b j+1 . . . b t a i+ i ...a k . (4.18) 

The multicommutator of (14.161) can thus be expanded in contraction schemes for 
operators at different times. An operator at time t\ contracts necessarily with 
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Figure 2: Illustration for the numbers r m , 
tractions. We see that they define a tree. 



, r 2 that occur in the choice of con- 



an operator at time t = 0; an operator at time t 2 contracts with an operator at 
time t r2 with r 2 = 0,1; . . . ; an operator at time t m contracts with an operator 
at time t Tm with r m = 0, . . . , m — 1. See Fig. El for an illustration. To a set 
of contractions corresponds a monomial of creation and annihilation operators, 
multiplied by anticommutators of contracted operators. 

The monomial of creation and annihilation operators is bounded in operator 
norm by 1. Contracted operators are estimated using (14. 15b . This yields a factor 
involving times, namely 

m 
3=1 

Second, one obtains a factor involving indices of Hermite functions for the con- 
tracted operators. An upper bound on this factor is obtained by writing a product 
over all indices, namely 



nnn(<&*+i)*(^*+i) 4 - 

j=0 k=l 1=1 

Here, </ - fc , q'- k E N d are indices for Hermite functions determined by the k-th 
element of Qj . 

It remains to estimate the number of contraction schemes, given r 2 , . . . , r m . 
We define 

ej = \{k : r k = j}\ + 1 - 5 jQ , 0<j<m. (4.19) 

Notice that 1 < e < m and 1 < ej < m — j ; + 1 if j ^ 0. ej is the number of 
operators at time tj that belong to a contraction and it is necessarily smaller than 
2Nj. Since there are 2Nj operators at time tj, the number of possible contractions 
is 

(2N,y. 



n 

j=0 



(2N 3 
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The above estimates could be improved by observing that many contraction schemes 
yield zero; namely, in the case where both operators are creation or annihilation 
operators; or if the spins are different; or if they belong to different reservoirs. It 
is not easy to take advantage of this, however. 

We now gather the above estimates to obtain the bound 

|[w(0,An-i(^]P< 

m 

x; i^ (Qr°\^ o) )ini^wr } ^; W) )i 

Qo,-,Q m j =1 

Nj d 

X 

7=0fc="l7=l ^oo>t ra >->t 1 >0 
m-l m-2 1 m m r / oArl/osr \l 

x e e -EndA^rn^y- <«<» 



nnrife+ i )"fe +i ) 1 / 

~_n J 00>t m > - >U >0 



=0r m _i=0 r 2 =0j=l j=0 



e 



3 1 



A sequence of numbers r 2 , . . . , r m can be represented by a graph with set of 
vertices {0, . . . , m}, and an edge between i and j whenever rj = i. This graph is 
a tree: there are m edges, and each vertex j ^ is directly connected to a vertex 
i < j, hence each vertex is eventually connected to 0. The numbers ej defined 
in (14.191) are then the incidence numbers of the tree — e 3 - is the number of edges 
containing the vertex j. This is illustrated in Fig. El We can symmetrize the bound 
by summing over all trees T with m + 1 vertices; this step will allow to deal with 
the time integrals. Reorganizing, we obtain 

\\[w(t),D m ^(t)a]\\dt< 

N d 

e \*nM n °\qo {no) )\ nn^ + + 1)' 

Q„ fe=l i=l 

/ *, . . . *. v n (i a «^) i/2 ^4^^ 

L>u>~>hx, T^mr (2iVo -<*>)! 



X 



^ ft I \ - , . ( ,\ i v.). 

j'=:i v Q, 

(4.21) 



j=l ^ Q, fc=l i=l ^ J 3 ' ' 



Let us focus on the time integrals. The integrand is a symmetric function of 
t\, . . . , t m , because of the sum over arbitrary trees. We can therefore extract a 
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factor 1/m!, at the cost of integrating over all positive times ti,...,t m without 
the ordering condition. Since there is no integral over t = 0, we have 

/oo poo _ poo m 

*,.../ dt m n (lAp^r^ / (lAff)" /2 di , (4.22) 

for any tree T. The last integral is equal to 87rd/ (d — 2). 

The number of trees with m + 1 vertices and incidence numbers e , . . . , e m is 
equal to 

/ m - 1 \ _ (m - 1)! 

Uo - 1, ei - 1, . . . , e m - 1 J ~ (e -l)!...(e m -l)! ; 



see for instance [Ber], Theoreme 2 p. 86. We sum over incidence numbers, using 

/e=0V e 



Ee=o( 2 e Ar )=^,andweget 



|| [W(f), ^i(t)a] ||* < ^ ( j^f E W^la^CQ^, Q; 



(iVo) ^/(Woh 




Qo 

jVo d 

xIIII(te + 1 )'^> + 1 )' 

fe=l i=l 

x E^i^ w ^' w )inn(^+ i ) | (^+ 1 ) 1 • (4 - 23) 

^ Q fc=l i=l J 

The last step consists in removing the Hermite functions. We fix TV, s^, 
s''^, r^ N \ r'( N \ and perform the summation over q( N ) and g'^. Using Cauchy- 
Schwarz, we obtain 

N d 

e |^(o (jv) .o /(jv) )|nn(9«+ i )'(^+ i ) i 

q (N)^ q ,(N) k=l i=l 

/ 2 N d x 1/2 

< e |^w (iv) ^ /(7v) )| nri(fe+ i ) 3 (^+ 1 ) 3 

fc=l i=l ' 

x(E^) ■<«<> 

The last factor on the right side is bounded by 2 dN . The first factor on the right 
side can be viewed as an expectation value of a certain operator expressed in the 
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basis of Hermite functions. Rewriting it in the a;-space representation, we find that 
it is given by the square root of the following expression 



2&dN 



R dN J R dN 



k=l i=l 



This motivates the use of the norm (14.41) and concludes the proof of Theorem 
I4~T1 □ 



We end this section by remarking that an estimate can be obtained that is 
invariant under space translations. Such an estimate follows by repeating the steps 
above with translates of Hermite functions. Recall that was defined in (14.101) 
by integrating 2dN Hermite functions centered at the origin. We can choose z E 
M. 2dN and translate the j-th function by Zj. Lemma H~3l still holds with translates 
of Hermite functions, so that the proof goes through without a change until (14. 241) . 
Since a Hermite function translated by z E R satisfies the differential equation 
(14.61) with (x — z) 2 instead of x 2 , one gets a bound where the differential operators 

in the norm (14.41) are translated by z E IR M . This holds for all z; let us introduce 

l " 

\m 



i" by 



// 

M 



r 

¥ U J 2 y m« 



/ dx^f(d^)f[(~ + (x h - y k ) 2 + lYf(x 
J r m dx k J 



(4.25) 



n ^ W 11/ 

I \M- 



This object is translation invariant but it is not a norm. We have 1 1 • \ \ M < 
Theorem l4. ll holds when ||A||' and ||W||' are replaced by \\A\\" and ||W||", whose 
definition is like (14.31) with || ■ ||^ instead of || • 



5 Explicit perturbative calculation of particle and 
energy currents 

In this section, we consider two reservoirs of non-relativistic non-interacting free 
spinless fermions. Such systems are a special case of the ones introduced in Sec- 
tion 

For explicit calculations, it is convenient to represent the system in Fourier 
space, see Subsection 15.11 since the one-particle energy operator t = —A is di- 
agonal in this representation. Subsection 15.21 is devoted to the calculation of the 
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particle and energy currents for tunnelling junctions, in the lowest non-vanishing 
order in W. This establishes the relation between the particle current and chemical 
potentials of the reservoirs, the current voltage characteristics. If the difference 
of chemical potentials A/x = // — fx 11 is small, then the particle current is pro- 
portional to the voltage drop A/i. This linear relation is known as Ohm's law. We 
calculate the (inverse of the) proportionality factor, which is called the resistance 
of the junction. Moreover, we explicitly verify that the entropy production rate is 
strictly positive, provided the two reservoirs are at either different temperatures or 
chemical potentials. 

Let us recall that the single particle Hilbert space (in the thermodynamic limit) 
is h — L 2 (R d ,dx), with x e IR d , d > 3 (see equation 93.21 ). The dynamics 
is determined by t = —A, see (13 ,3b . For each reservoir, we take the particle 
number to be the only conservation law. Recall that for tunnelling junctions, the 
interaction W commutes with the total particle number operator, N <g> 11 + 11 ® N, 
while for thermal junctions, W commutes separately with N <g> 1 and 11 ® N; see 
equations (13.411) and (13.401) . 

To quantify the interaction, we introduce two coupling constants, g and £, and 

set 

oo 

W = g^ N W N . (5.1) 

N=l 

Let J7fc,z, Vk,i denote the term of order g k £ l of the particle current J (see (12.861) . 
(I2~78t dZHB ) and the energy current V (see (l2~85l) . (127771 . (l2~82T ). Accordingly, 
we define £k,i, where £ is the entropy production rate in (12.691) and (12.75b . We 
give now explicit expressions for some lower order terms of the currents and the 
entropy production rate. The calculations are presented in Subsection 5.2. 

Tunnelling junctions. The lowest order terms of the particle current are given 

by 

Ji,i = Ji,2 = 0, (5.2) 
J 2i2 = 2vr / dkdUilkf-ll^lw^i-kJ^^lJVfipuity-pjik)), 

JR 2d 

(5.3) 

where k,l E M. d , w is the Fourier transform of w, and the function p r (k) is defined 
as 



FMU February 7, 2008 



43 



We obtain for the energy current the expressions 

V 1A = Pi, 2 = 0, (5.4) 
V 2>2 = 2n [ rffc^|fc| 2 5(|fc| 2 -|Z| 2 )|«; 1 ((-fc,//),(Z,/))| 2 (p // (fc)-P/(fc))- 

jR 2d 

(5.5) 

Assuming that w 1 is not identically zero, the above formulas show the following 
qualitative behaviour of the flows. 

- If (P^p 1 ) = (P n ,p n ) then J 2>2 = V 2 , 2 = 0. The flows vanish if both 
reservoirs are at the same temperature and chemical potential. 

- If p 1 = (j, 11 and (3 1 > /3 H then pn(k) - p^k) > 0, for all k. Consequently, 
1/2,2, ^2,2 > 0. At constant chemical potential, there is a particle and energy 
flow from the hotter to the colder reservoir. 

- If p 1 > p n and p 1 = p n = p, then p n {k) - Pl {k) < 0, for all k. 
Consequently, J 2 ,2, ^2,2 < 0. At constant temperature, there is a particle 
and energy flow from the reservoir with higher chemical potential to the 
reservoir with lower chemical potential. 

Ohm's law and the resistance of the junction. Suppose that (3 1 = [3 11 and 
p 1 = p, p 11 = p + Ap, with Ap small. Retaining only the leading order in Ap 
in the expression of the particle flow yields 



Ap 

rJpTPY 



J2,2 w 777— —777, (5.6) 



where the resistance R(p, (3) is determined by 

R(p,P) =2n? J^dkdl5{\k\ -\l\ ) (e , ( | fc | 2 - M) + 1)2 • 

(5.7) 

We refer to Subsection 5.2 for a qualitative discussion of the resistance, in three 
dimensions, d = 3. 



Onsager reciprocity relations. Let us study the interdependence of the flows 
near equilibrium. The relevant parameters are the difference of the inverse temper- 
atures, and the difference of the chemical potentials divided by the temperature. 
Precisely, we set (3 1 = (3; f3 H = (3 - Ap; v = p T p J ; Av = (3 1 p 1 - p H p H . We 
consider the flows to depend on P, v, A/3, and Av. 
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One easily checks that 





dAp 
d 



[pn{k] - p/(fe 
[p n (k) - pi{k 



a/3=a,=o ( e «l fc l 2 ^) + l) 2 

& 



k\ 



o/3(|fc| 2 -M) 



A/3=A„=0 ( e «|fe| 2 -M) + l)2- 



(5.8) 



The first partial derivative is taken at constant (3, v, and Av; the second partial 
derivative is at constant (3, v, and Af3. Then from (15.31) and (15.51) we observe that 



2.2 



<9Az/ 



&7s 



2.2 



A/3=Ai/=0 



dA(3 



A/3=Av=0 



(5.9) 



This is an Onsager reciprocity relation and we see that it holds at lowest order. 

Entropy production rate. Recall that V H = —V 1 (equation (12.771) ) and J 11 ■ 
—J 1 (equation (I2.78!) ). hence 

S = (f3 I - (3 n )V - (/?V - /3 n p n )J. 

Using the above expressions for J k i and Vk ,i, we obtain 

£1,1 = &i,2 — 0? 

^2,2 = -f3 II )V 2 , 2 -{(3 I p I -I3 n p n )j %2 

= 2vr I dkdl5(\k\ 2 - \l\ 2 )\wi{(-k, II), (I, I))\ 2 



x 



( e ^(|fc| a -^) + i)( e /3"(|fc| 2 -M") + i) 



The numerator of the fraction is of the form (x 1 — x IT )(e x — e xI1 ), with x r 
(3 r (\k\ 2 — p r ), hence it is strictly positive unless x 1 = x 11 . We assume that 



/ 



dlSdkf-ll^lwtd-kJlKl,!))? 



(5.10) 



does not vanish for all k E M. d . Then £ 2j2 is strictly positive unless (3 1 (\k\ 2 — p 1 ) = 
(5 n {\k\ 2 — p 11 ) for all k in the support of (15.101) . This shows that £ 2<2 is strictly 
positive unless (/3 7 , p 1 ) = {(3 11 , p 11 ), in which case £ 2i2 vanishes. 



Thermal junctions. The particle current is zero, a thermal junction allows 
only for an exchange of heat between the two reservoirs. Since Wi = 0, the lowest 
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order term which is nonvanishing is "P 2j4 . Without loss of generality, we take the 
coupling function w 2 to be of the form 

W 2 ((fcl, n), (fc 2 , ^2); (Zl, (Z 2 , S 2 )) = $ri,I $r 2 ,n 8 si ,i S S2>n w 2 {k 1 , k 2 , h, Z 2 ). 
A somewhat lengthy but straightforward calculation yields 

P 2>4 = 2n J dk x dk 2 dhdl 2 (Klfcil 2 + IN 2 " l*i| 2 " N*) 

x|^ 2 (-fc 1 ,-fc 2 ,Z 1 ,Z 2 )| 2 (|fc 1 | 2 -|Z 1 | 2 )p / (Z 1 ) P// (^)(l-p/(fci)-p//(fc 2 )), 

(5.11) 

from which we obtain the following qualitative discussion. 

- If (/3 1 , //) = (f3 n , p n ), then pj = p n , and by switching Zi <-> Z 2 , fci fc 2 
in the integral, we see that "P 2 4 = 0. 

- By splitting the integral in (15.1 II) into a sum of two integrals over the regions 
x(|fci| 2 > |Zi| 2 ) and x(l&i| 2 < |Zi| 2 ), and switching ki <-> Zi, k 2 <-> Z 2 , we 
can rewrite 



p 24 = 2tt y dk t dk 2 dl x dl 2 8{\k^ + |fc 2 | 2 - |Zx| 2 - |Z 2 | 2 ) 

x|«; 2 (-fc 1 ,-fc 2 ,Z 1 ,Z 2 )| 2 (Ifcxl 2 - |Zx| 2 ) X (|fci| 2 > 
x{p / (Zi)p // (Z 2 )[l - pi{kt) - pii{k 2 )} 
-pi(h)pn(l 2 )[l - pj(fei) - Pn(k 2 )]} 

= 2ix [ dk x dk 2 dhdl 2 8{\k x \ 2 + |fc 2 | 2 - \h\ 2 - |Z 2 | 2 ) 



x\w 2 (-k u -k 2 ,l l7 l 2 )\ 2 (Ifcxl 2 - |Zx| 2 ) X (\ki\ 2 > \h\ 2 ) 
x{p 7I (Z 2 )[l - pn(k 2 )][pi(li) -pi(fei)] 
-p/(*i)[l-p/(Ii)]M*a)-P//('2)]}. 

The first product in the round brackets { } is strictly positive and tends 
to zero, as (3 H — > 00 (because in the limit /3 r — > 00, ^(fc) tends to the 
characteristic function x{ I k \ 2 < p T )). The second term in the round brackets 
is strictly negative and tends to zero, as (3 1 — > 00. We conclude that V 2 >4 < 
if (3 1 < 00, and /J 77 is large enough; as expected! 
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5.1 Fourier representation 

The creation and annihilation operators in the Fourier representation are defined 
by 

a\k,r) = (2n)- d l 2 [ dxe ihx a*{x,r), (5.12) 

JR d 

a{k,r) = (2ir)- d/2 [ dx e - tkx a(x, r), (5.13) 

J~R d 

where k G M d , r — I, II; compare with (13.161) and (13 .22ft . The dynamics of 
a* (k, r) and a{k, r) is given by 

a r t (a*(k,r)) = e itJt a*(k,r), a 7 t {a{k,r)) = e~ lujt a{k,r), (5.14) 

where 

LU = Uj(k) = \k\ 2 . 

The operators H r , N r , W N defined in (13 . 19b . (13.20b . (13.341) are represented in 
Fourier space (and in the thermodynamic limit) by 

H r = I dku(k)a*(k,r)a(k,r), (5.15) 

JR d 

N r = I dka*{k,r)a{k,r), 

W N = JdK^^dL^(^)w N (-^,L^)a{L^), 

where we introduce notation analogous to (I3.24I) - (I3.31I) . For = (Ki, . . . , K N ), 
Kj = (kj, rj ) E R d x {/,//}, we set -K™ := {-K u . . . , -K N ), where 
— Kj = (—kj, Tj). The symbol ~ denotes the Fourier transform, i.e. 

= (2ir)- dN [ d Xl ---dy N e- i{klXl+ - +lNyN) w N (X^ N \Y^). 

We recall some properties of the state cu° defined in (12.601) . which is given by 

uj° = up <g> ujpu tfl n, (5.16) 

where upr )ft r is the equilibrium state of reservoir r — 1,11 in the thermodynamic 
limit; see also Theorem l3.ll The two-point function of LOpr^r is 

upr tfl r(a*{k,r)a(l,r)) = 5{k - l)p r (k), where p r (k) = 2 _ . 
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The average of a monomial in n creation and m annihilation operators is zero 
unless n = m, in which case it can be calculated recursively from the formula 

(n n 
n a *( fci ' r )ii a( ^' r ) 
i=\ j=i 

n / n n 

= ^2(-~L) n ~ p upr iA ,r(a" '(k u r)a(l Pi r)) upr^r I JJa*(fc;,r) JJ a(Z i ,r) 

p=l \i=2 j=l t j^ p 

For details, we refer to [BR]. We are now ready for explicit calculations of the 
currents. 

5.2 Calculations for tunnelling junctions 

Particle current and resistance. The particle current 

J = ^statHfiV <g> 1, W]) = uj°(a+(-i[N <g> 1, W])) (5.17) 
has been introduced in (12.861) . (12.781) . see also (12.811) . We set h = 1. We see from 



the Dyson series expansion of er + , see (14.11) . and the definition of Jk,i ( see a f ter 
(ED), that 

J 1A = -iu°([N ®l,Wi]), 
Ji, 2 = -io;°([iV®]l,W 2 ]), 

/■oo 

J 2i2 = / tftu; ([Wi(t)jW®ll > Wi]]). (5.18) 
It is not difficult to verify that 



N 



[N®1,W N ] = I dK^dL^J2{6 rj>I -6 r/j>I ) 

J 3=1 

x a *(K^)w N (-K( N \ LW)a(LW), (5.19) 

where 5 is the Kronecker symbol and = (Li, . . . , L/v)> £j = r'X Using 
that a; is invariant under A i-» e lsN ' Ae~ lsNl , we find that 

w°([JV®ll,Wi])=0, for all/. (5.20) 

Thus = Ji j2 = 0. Next, we calculate J^- Recalling that Wi{t) = a°_ t (W) 
and equation (15.141) . we write 

\Wi(t),[N®l,Wi]]= V {S s ,i-6 sl> i) I dk dl dk' dV e-^l 2 '!'! 2 )' 

r,r',s, S '=I,II ' Jmid 
xwtd-k, r), (/, r')) wxd-k', s), (V, s')) [a*(k, r)a(l, r'), a*{k', s)a(l', s% 
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We expand the commutator on the right side and apply the state u° to obtain 

poo p 

J 2 , 2 = dt dkdlw 1 ({-k i II),{l > I))w 1 ({-l > r),(k,II)) 

JO JR 2d 

x | e -i(| fe P-|^) t + eWHO } ( pn (k) - pj(l)) . 

Because W\ is selfadjoint, we have the relation 



Gi((-k, II), (I, I)) = wx((-l, I), (k, II)); 

using this and the formula dt e 1Tt = 2n5(r), one sees that (15.31) holds. It is 
useful to keep in mind that 



Pii(k) - Pl (k) 



e /3^(|fc| 2 -^) _ e /3"(|fc| 2 -^") 



3 ^(|fe|2-^) 



+ 1) ( e ^(|fcP-M») + i) 



The resistance. Let (3 1 = (3 n and p 1 = p, p 11 — p + Ap, with Ap small. We 
expand 

PlI (k) - pj(k) = (3 A/i e m?-»)^m?-») + !)-2 + 0((A/i) 2 ). 

Retaining only the first order in A/i in equation (15.31) gives (15.61) and (15.71) . Let 
T = 1/(3 denote the temperature and assume that p > 0. We see that R(p, (3) ~ 
T, as T — > oo. Next, we examine the dependence of the resistance on T, for small 
T, in three dimensions, and where w\ is a radial function in both variables (i.e. W\ 
depends only on |fe| and \l\). We then have 

poo pP{r-ii) 

Rip^r^Sn'P / drr|^ 1 (( v ^,//),( v ^,/))| 2 



(e/3(r-M) + l) 2 ' 



The fraction in the integral equals — (3 1 d r (e^ r w + 1) 1 and it follows that 

f°° 1 
R(p,(3)-^J o drf'(r) e ^ ) + v (5.21) 

where /' denotes the derivative of the function 

/(r) :=87r 3 r|e i (( v ^,J/),( v ^,/))| 2 . 
Let us split the integral in (15.211) as 



R (p,pr = nr)-^-——+ f(r) 



Pfl ^ POO 



1 

+ : 

(5.22) 
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Apply the change of variables t = —f3{r — fi) and t — /3(r — fi) in the first and 
second integral on the right side of (15.221) . repectively. Then one has 



R(ji,P)- l = f(ji) 
1 f°° 

+ J dt {fit IP + n) - f(-t//3 + n) X (t < (e* + I)" 1 
and using the mean value theorem, 

R^vr 1 = m + j 2 QHcft r^)^ + o( e -^)) , 

for some & G [—t/P + fi, t/ (3 + /x] and where the exponentially small remainder 
term comes from removing the cutoff function x(t < (3jj). Retaining the main 
term {J3 — > oo) yields 

'2 



and consequently, 



At zero temperature, the resistance has the value R(fi, oo) = and it in- 

creases or decreases with increasing T according to whether /"(/x) < Oor/"(/x) > 
0. 

Energy current. The energy current 

V = 0J 8ta t(-i[H ®1,W]) = u°(a+(-i[H ® 1, W])) 

has been introduced in equations (I2.85I) . (12.77b . see also (I2.82I) . We set /t, = 1. 
Using the CAR and expression (15.151) for H, one obtains 

r N 

[H ® 11, W^v] = / dK {N) dL w ^(|fejf«y PiJ j - |Zi| 2 ^.,/) 
^ j=i 

x a* ( ATW )«W ( -^ (7V) , L W ) a(L W ) , 

and it is readily verified that V\ \ = V\ i = 0, and a similar calculation as for the 
particle current shows that 

POO 

P 2i2 = / dt [H ® 1, W x ]]) (5.23) 

Jo 

is given by (I5.5I) . 



FMU February 7, 2008 



50 



References 

[Af] J. Affleck: Fields, strings and critical phenomena, (Les Houches 1988), 
Amsterdam North Holland, 1990 

[ACF] A. Yu. Alekseev, V.V Cheianov, J. Frohlich: Universality of Transport 
Properties in Equilibrium, the Goldstone Theorem, and Chiral Anomaly, 
Phys. Rev. Letters 81 (1998), 3503-3506 

[AP] W. Aschbacher, C.-A. Pillet: Non-Equilibrium Steady States of the XY 
Chain, preprint, mp-arc, 2002 

[ASS] J. E. Avron, R. Seiler, B. Simon: Charge deficiency, charge transport and 
comparison of dimensions, Comm. Math. Phys. 159 (1994), no. 2, 399-422 

[Bee] C. W. J. Beenakker: Random-matrix theory of quantum transport, 
larXiv:cond-mat/9612179k l. 19 Dec. 1996 

[Ber] C. Berge, Principes de Combinatoire, Dunod, Paris, 1968 

[BES] J. Bellissard, A. van Elst, H. Schulz-Baldes: The noncommutative geom- 
etry of the quantum Hall effect. Topology and physics., J. Math. Phys. 35 
(1994), no. 10, 5373-5451 

[BFS] V. Bach, J. Frohlich, I. M. Sigal: Return to Equilibrium, J. Math. Phys. 41 
(2000), No. 6, 3985-4060 

[BLR] F. Bonetto, J.L. Lebowitz, L. Rey-Bellet: Fourier's Law: A Challenge to 
Theorists, Mathematical physics 2000, 128-150, Imp. Coll. Press, London, 
2000 

[BM] D. D. Botvich, V. A. Malishev: Unitary equivalence of temperature dy- 
namics for ideal and locally perturbed Fermi-gas, Comm. Math. Phys. 91 
(1983), no. 3, 301-312 

[BR] O. Bratteli, D. Robinson: Operator Algebras and Quantum Statistical Me- 
chanics I, II, Texts and Monographs in Physics, Springer Verlag, 1987, 1997 

[Ca] H. B. Callen: Thermodynamics and introduction to thermostatistics, John 
Wiley & Sons, second edition, 1985 

[CFKS] H.L. Cycon, R.G. Froese, W. Kirsch, B. Simon: Schrddinger Operators, 
Texts and Monographs in Physics, Springer Verlag, 1987 

[CNP] C. Connaughton, A. C. Newell, Y. Pomeau: N on- stationary spectra of 
local wave turbulence, preprint (2003) 



FMU February 7, 2008 



51 



[CPW] A. M. Chang, L. N. Pfeiffer, K. W. West: Observation of Chiral Luttinger 
Behaviour in Electron Tunneling into Fractional Quantum Hall Edges, Phys. 
Rev. Letters, 77 (1996), no. 12, 2538-2541 

[CWCPW] A. M. Chang, M. K. Wu, C. C. Chi, L. N. Pfeiffer, K. W. West: 
Plateau Behaviour in the Chiral Luttinger Liquid Exponent, Phys. Rev. Let- 
ters, 86 (2001), no. 1, 143-146 

[DFG] Dirren, S., ETH diploma thesis winter 1998/99, chapter 5; (written under 
the supervision of J. Frohlich and G.M. Graf) 

[EPR] J. -P. Eckmann, C.-A. Pillet, L. Rey-Bellet: Entropy production in nonlin- 
ear, thermally driven Hamiltonian systems, J. Stat. Phys. 95 (1999), no. 1-2, 
305-331 

[FMSU] J. Frohlich, M. Merkli, S. Schwarz, D. Ueltschi: Statistical Mechanics 
of Thermodynamic Processes, submitted to World Scientific Press 

[FP1] J. Frohlich, B. Pedrini: New applications of the chiral anomaly, Mathemat- 
ical physics 2000, 9-47, Imp. Coll. Press, London, 2000. 

[FP2] J. Frohlich, B. Pedrini: Axions, quantum mechanical pumping, and 
primeval magnetic fields, in A. Cappelli and G Mussardo (eds.), Statistical 
Field Theories, 289-303, Kluwer Academic Publishers, 2002 

[FST] J. Frohlich, U. M. Studer, E. Thiran: Quantum theory of large systems 
of non-relativistic matter, Geometries fluctuantes en mecanique statistique 
et en theorie des champs (Les Houches, 1994), 771-912, North-Holland, 
Amsterdam, 1996. 

[Ha] R. Haag: Local Quantum Physics, Texts and Monographs in Physics, 
Springer Verlag, 1992 

[He] Hepp, K.: Rigorous results on the s — d model of the Kondo effect. Solid 
State Communications, 8 (1970), 2087-2090 

[JP1] V. Jaksic, C.-A. Pillet: On entropy production in quantum statistical me- 
chanics, Comm. Math. Phys. 217 (2001), no. 2, 285-293 

[JP2] V. Jaksic, C.-A. Pillet: Non-equilibrium steady states of finite quantum stys- 
tems coupled to thermal reservoirs, Comm. Math. Phys., 226 (2002), No. 1, 
131-162 

[La] R. B. Laughlin: Quantized Hall Conductivity in Two Dimensions, Phys. Rev. 
623(1981) 5632-5633 



FMU February 7, 2008 



52 



[LS] L. S. Levitov, A. V. Shytov: Tunneling in a bosonized Fermi liquid, 
larXiv:cond-mat /95 10006 . 1995, release of September 2002 

[LSH] L. S. Levitov, A. V. Shytov, B. I. Halperin: Effective action of a compress- 
ible QH state edge: application to tunneling, arXiv:cond-mat/0005016^2, 
July 2001 

[Ma] G. D. Mahan: Many Particle Physics, Plenum Press, New York, 1981 

[MO] T. Matsui, Y. Ogata: Variational principle for non-equilibrium steady 
states of the XX model, preprint (2003) 

[Rob] D. Robinson: Return to Equilibrium, Comm. Math. Phys. 31 (1973), 171— 
189 

[RSII, III] Reed, M., Simon, B.: Methods of Modern Mathematical Physics, Vol. 
II (Fourier Analysis, Self-Adjointness), III (Scattering Theory), Academic 
Press, New York 1975, 1979 

[Rul] Ruelle, D.: Natural nonequilibrium states in quantum statistical mechan- 
ics, J. Stat. Phys. 98 (2000), no. 1-2, 57-75 

[Ru2] Ruelle, D.: Entropy production in quantum spin systems, Comm. Math. 
Phys. 224 (2001), no. 1,3-16. 

[Ru3] Ruelle, D.: Statistical Mechanics. Rigorous results, Reprint of the 1989 
edition. World Scientific Publishing Co., Inc., River Edge, NJ; Imperial Col- 
lege Press, London, 1999 

[SGJE] L. Saminadayar, D. C. Glattli, Y. Jin, B. Etienne: Observation of the e/3 
Fractionally Charged Laughlin Quasiparticle, Phys. Rev. Lett. 79 (1997), 
2526-2529 

[SRI] M. Sigrist, T. M. Rice: Paramagnetic Effect in High T c Superconductors 
-A Hint for d-Wave Superconductivity, Journal of The Physical Society of 
Japan, 61 (1992), no. 12, December, 4283-4286 

[SR2] M. Sigrist, T. M. Rice: Unusual paramagnetic phenomena in granular 
high-temperature superconductors - A consequence of d-wave pairing?, 
Rev. Mod. Phys., 67 (1995), no.2, April, 503-513 

[TKNN] D. J. Thouless, M. Kohmoto, M. P. Nightingale, M. den Nijs: Quantized 
Hall Conductance in a Two-Dimensional Periodic Potential, Phys. Rev. Lett. 
49 (1982), 405-408 



FMU February 7, 2008 



53 



[vW] B. J. van Wees, H. van Houten, C. W. J. Beenakker, J. G. Williamson, L. P. 
Kouwenhoven, D. van der Marel: Quantized Conductance of Point Contacts 
in a two-dimensional Electron Gas, Phys. Rev. Lett. 60 (1988), 848-850 



